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ABSTRACT 


Methods  were  developed  to  predict  the  boundary  layer  and  wake  characteris¬ 
tics  of  symmetric  uodies  rrv  •  •/  ’  at  supersonic  speeds.  The  methods  devel¬ 
oped  for  predicting  boundary  lr.yer  characteriit  .cs  are  applicable  to  both 
laminar  and  turbulent  flow  over  axisymmetric  or  two  dimensional  bodies  and 
can  be  used  as  an  approximation  for  the  flow  over  quasiaymmetric  bodies. 
The  wake  predictive  methods  were  developed  for  axisymmetric  flow  and  are 
applicable  to  both  laminar  and  turbulent  wakes.  The  aforementioned  methods 
were  combined  and  programmed  such  that  wake  profiles  may  be  determined 
based  only  on  inputs  of  forebody  geometry  and  free  stream  conditions.  Also 
presented  is  a  method  for  calculating  the  drag  of  0Ar.ASONIC^  parachutes 
using  the  flow  field  immediately  ahead  of  the  parachute  as  the  free  stream 
conditions.  An  analysis  of  the  effects  of  various  geometric  and  free  stream 
parameters  on  PARASONIC  parachute  performance  wag  conducted  and  trends 
determined. 

The  distribution  of  this  abstract  is  unlimited. 
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U'/U  =  free  stream  turbulence 
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T  -  temperature 


AFFDL-TR-67-192 
Volume  II 


u  =  x-direction  velocity  component 
Vjjj  =  free  stream  velocity 

v  =  radial  direction  velocity  component 
X  =  downstream  distance  from  forebody  base 
=  rear  stagnation  point  location 

Y  =  vertical  distance  from  forebody  centerline  (+  above 
and  -  below) 

Z  *  horizontal  distance  from  forebody  centerline  {+  right 
and  -  left,  looking  upstream) 

y  =  ratio  of  specinc  heat  capacities 

A  ~  transformed  wake  radius  or  thickness 

S  =  wake  radius  or  thickness 

ip  =  wake  diameter,  6^  =  26 

6*  ~  wake  displacement  thickness 

rj  -  nondimensional  transformed  radial  coordinate 
0  a  momentum  thickness 

0-j.  =  enthalpy  thickness 
u  =  viscosity 

£  =  wake  form  parameter  in  Equation  65 
fi  =  gas  density 

Subscripts 

B  =  forebody  base  where  boundary  layer  separate* 

D  =  (orehody  cLamcter 
L  =  locai  conditions 
t  =  total  condition 
s  =  forebodv  surface  condition 
TR  =  transition  condition 


AFFDL-TR-67-192 
Volume  II 


w  =  wall  condition 

x  =  distance  from  rear  stagnation  point 

oo  =  freestream  condition 

<5—  =  wake  diameter  at  transition 
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DRAG  PREDICTION  {SECTION  V);  PARASONIC  PARACHUTE  PERFORM¬ 
ANCE  (SECTION  IV);  WIND-TUNNEL  TESTING  AND  DATA  REDUCTION 
{APPENDIX  V) 
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Subsc  ripts 


P  =  pressure 
q  =  dynamic  pressure 
R  -  maximum  radius  of  parachute 
Re  =  Reynolds  number 
r  a  radius  of  ring  of  differential  width  ds 
r^  ~  inlet  radius  of  ring  of  constant  pressure  differential 

r  =  exit  radius  of  ring  of  constant  pressure  differential 


o 
WR 
X 
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Z  = 


a  = 


3  = 
x  = 
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width  of  radials 

downstream  distance  from  forebody  base 

vertical  distance  from  forebody  centerline  (+  above 
and  *  below 

horizontal  distance  from  forebody  centerline  (+  right 
and  -  left,  looking  upstream) 

angle  between  chord  and  skirt  axis 
«/  M2  -  1 

geometry  porosity 


=  total  porosity 

-  angle  between  a  line  tangent  to  parachute  canopy  and 
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b  = 


i  = 


L  - 


torebody  base  or  breathing 
constructed  or  ,J nopy 
mesh  opening 
drag  oscillation 
inlet 

local  conditions 
stagnation  conditions 
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p  =  projected 
t  =  total  conditions 
exp  =  experimental 
ext  =  external 
int  =  internal 

a  =  canopy  oscillation  angle 

1  =  conditions  forward  of  Bhock  wave 

2  =  conditions  aft  of  shock  wave 
aa  =  free  stream  conditions 

FOREBODE  FLOW  FIELD  PROGRAM  (APPENDIX  II) 
i  =  axial  indexing  parameter 
M  =  Mach  number 
p  =  pressure 

/ 

R  =  average  base  radius 

r  -  radial  coordinate  normal  to  body  centerline 
y  =  ratio  of  specific  heats 

0  =  angle  of  the  forebody  surface  with  the  body  axis  in 
each  f6  plane 

0  =  critical  cone  angle  for  attached  shock 

c 

^  =  Reference  plane  passing  through  body  centerline 

Subscripts 

e  =  Conditions  at  boundary  layer  edge 
o  =  local  stagnation  conditions 
s  =  nose  stagnation  conditions 
co  =  free  stream  conditions 


(Reverse  is  blank) 
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SECTION  I 
INTRODUCTION 


As  a  portion  of  the  Establishment  of  an  Unsymmetrical  Wake  Test  Capability 
tor  Aerodynamic  Decelerators  (EUREKA)  program  sponsored  by  the  Air 
Force  Flight  Dynamics  Laboratory  (AFFDL),  a  study  was  performed  to  ana¬ 
lyze  performance  characteristics  of  supersonic  parachutes  operating  in  the 
wakes  of  both  axisymmetric  and  asymmetric  forebodies  and  to  develop  the 
capability  to  analytically  predict  the  performance  of  a  parachute  operating 
in  the  wake  of  an  axisymmetric  forebody.  The  results  of  that  study  are  pre¬ 
sented  in  this  volume, 
f 

In  order  to  predict  the  parachute  performance”,  the  properties  of  the  wake 
flow  field  in  which  the  parachutes  operate  and  the  forebody  flow  field  proper¬ 
ties  that  establish  the  development  of  the  wake  flow  field  had  to  be  determined. 
The  three  basic  tasks  performed  were: 

1.  Analysis  of  the  flow  field  and  boundary  layer  character¬ 
istics  of  the  axisymmetric  forebodies 

2.  Determination  of  the  viscous  and  inviscid  wake  flow  field 
characteristics  behind  these  forebodies 

3.  Prediction  of  the  drag  performance  characteristics  of 
the  parachute  operating  in  the  wake  of  such  forebodies 

Presented  in  Section  II  is  a  method  for  determining  the  flow  properties  over 
the  forebody  and  the  boundc.  •  ’ayer  momentum  thickness  and  momentum  de¬ 
fect  on  the  body.  The  invia-iu  i!  >.v  properties  necessary  for  the  boundary 
layer  analysis  are  computed  using  Newtu.ian  and  tangent  cone  theories, 
where  the  tangent  cone  theory  is  based  on  the  Taylor-Maccoll  conical  flow 
solution.  The  laminar  compressible  boundary  layer  characteristics  are  de¬ 
termined  by  using  a  IocjI  similarity  solution,  and  the  turbulent  boundary 
layer  is  analyzed  by  using  an  integral  method.  An  approximate  criterion  for 
the  boundary  layer  transition  point  also  is  given. 

Presented  in  Section  III  is  a  momentum  integral  solution  for  laminar  and  turbulent 
viscous  wake  characteristics  and  a  mass  flow  balance  solution  for  the  invis¬ 
cid  wake  characteristics.  The  viscous  wake  theory  is  based  on  the  assump¬ 
tion  that  the  details  of  the  base  flow  and  free  shear  layer  regions  may  be 
largely  ignored,  thus  allowing  the  wake  momentum  defect  to  be  equated  to  the 
forebod)  momentum  defect.  Boundary  layer  type  equations  then  can  be  applied 
to  the  wake  flow  region  downstream  of  the  wake  rear  stagnation  point.  Also 
presentee  in  Section  111  is  a  method  for  determining  inviscid  wake  character¬ 
istics  by  a  mass  flow  balance  that  is  performed  with  the  forebody  bow  shock 
shape. 

Based  on  the  viscous  and  inviscid  wake  solutions  along  with  the  computed 
forebody  characteristics,  comparison  cf  predicted  wake  profiles  has  been 
made  with  evperimental  wake  profile  data  obtained  from  EUREKA  wind- 
tunnel  tests  behind  an  Arapaho  C  with  nose  cone  forebody.  These  correlations 
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indicate  chat  excellent  agreement  between  wake  theory  and  experiment  can  be 
obtained  using  the  methods  presented  in  this  volume  The  wake  methods  of 
solution,  as  well  as  the  forebody  methods  of  solution,  have  been  programmed 
in  FORTRAN  IV  for  use  on  the  IBM  360,  Model  40,  and  IBM  7094  computers. 
This  computer  program  was  used  in  making  ail  the  wake  comparisons  shown. 


Section  IV  contains  an  analysis  of  the  effects  of  Mach  number,  Reynolds  num¬ 
ber,  forebody  shape,  canopy  location,  and  canopy  size  on  PARASOMIC*  para¬ 
chute  performance. 


Section  V  presents  a  method  for  calculating  the  pressure  drag  coefficient  of 
a  PARASON1C  parachute  operating  in  the  wake  of  an  axisymmetric  forebody. 
Comprrisons  of  calculated  and  experimentally  determined  drag  coefficients 
are  presented  also  and  show  good  agreement. 

Other  sections  of  this  volume  contain  conclusions;  recommendations  for 
future  work,  and  appendixes  giving  the  details  of  the  flow  field  analyses  and 
w  ad  tunnel  tests. 


This  report  does  not  cover  all  primary  program  objectives.  Volume  I  pre¬ 
sents  analyses  and  preliminary  design  efforts  to  determine  feasibility  of  re¬ 
configuring  a  free  flight  decelerator  test  vehicle  with  inflatable  appendage  in 
order  to  simulate  the  w^.^.  of  a  lifting  body. 

Volume  III  contains  a  tabulation  of  wake  survey  and  body  suria  e  pressure 
cata  obtained  from  wind  tunnel  tests  conducted  during  the  EUREKA  program. 
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SECTION  II 


BOUNDARY  LAYER  ANALYSIS 


BOUNDARY  LAYER  CHARACTERISTICS 


In  investigating  the  various  characteristics  of  the  viscous  wake  trail¬ 
ing  a  typical  aerodynamic  body,  one  of  the  dominant  features  is  integrated 
momentum  defect.  The  wake  momentum  defect,  hi  ,  can  be  expressed  as 


n 


M 


2n 


pu(u.  -  u)  rd^  dr 


and  for  a  wake  in  which  there  is  little  or  no  static  pressure  gradient, 
the  momentum  defect  remains  invariant  as  the  wake  develops  in  the 
downstream  direction.  Furthermore,  by  neglecting  the  details  of  the 
base  flow  r,d  shear  layer  regions  and  assuming  that  no  external  work 
is  done  on  .tie  wake  fluid  in  these  regions,  the  wake  momentum  defect, 

M  ,  may  be  evaluated  by  equating  it  to  the  integrated  boundary  h  it 
momentum  defect,  M,  at  the  trailing  edge  of  the  forebody.  The  equation 
at  momentum  defects  provides  a  direct  link  between  the  viscous  wake 
properties  and  the  forebody  configuration,  and  any  attempt  to  predict 
the  properties  of  the  viscous  wake  behind  a  body  must  begin  with  some 
knowledge  of  the  forebody  trailing  edge  boundary  layer  momentum  defect. 


This  section  describes  the  development  of  a  method  for  calculating  the 
growth  of  the  boundary  layer  momentum  thickness  and  momentum  defect 
over  an  aerodynamic  body  of  arbitrary  shape.  In  doing  this,  an  exten¬ 
sive  literature  survey  was  made  of  the  available  literature  on  compress¬ 
ible  boundary  layer  flow.  The  method  finally  chosen  for  use  on  the  pro¬ 
gram  is  representative  of  the  state  o i  the  art  of  present  boundary  layer 
theory.  The  resultant  procedures  developed  yield  solution'  for  two- 
dimensional  or  axisymmetric  flow  and  are  good  approximations  for  quasi- 
two -dimensional  or  quasi-symmetric  flows.  The  method  of  analysis  as 
developed  for  laminar  and  turbulent  flow  has  been  programmed  for  use 
on  high  .peed  digital  computers  and  has  been  utilized  to  determine  the 
momentum  defect  and  momentum  thickness  on  the  Arapaho  C  vehicle  with 
nose  cone  and  on  a  blunted  elliptical  cone  configuration. 


BOUNDARY  LAYER  CALCULATION  METHODS 


a.  Laminar  Boundary  Layer 


Kemp,  Rose,  and  Detra  (Reference  1)  have  studied  the  compressible 
laminar  boundary  layer,  including  effects  due  to  pressure  gradient  and 
due  to  chemical  dissociation.  The  former  effect  has  been  investigated 
through  the  study  of  similar  solutions  that  may  be  obtained  for  cases  in 
which  the  boundary  layer  pressure  gradient  parameter,  is  constant 
with  respect  to  surface  distance.  Although  for  most  cases  of  interest 
0  is  not  constant  over  the  entire  surface  of  a  body,  it  has  been  generally 
accepted  that  the  results  of  su.h  calculations  still  may  be  applied  in 
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engineering  calculations  if  the  local  value  of  the  pressure  gradient  para¬ 
meter,  j3,  is  used  (References  2  and  3).  This  application  is  called  the 
concept  of  "local  similarity.  "  Lees  in  Reference  4  has  shown  that  for 
highly  cooled  boundary  il-c  ;ffect  of  the  parameter  /3  is  small  and 

may  even  be  neglected  entirely  for  heat  transf'***  calculations.  However, 
for  skin  friction  calculations  and  for  cases  where  the  wall  temperature 
is  comparable  to  the  free  stream  total  temperature,  0  still  can  have  an 
important  effect.  Within  the  assumption  of  local  similarity,  however, 
this  influence  rnay  be  accounted  for  by  using  the  local  value  of  j3. 

As  was  indicated,  Kemp,  Rose,  and  Detra  also  have  included  effects  due 
to  dissociation  by  using  the  binary  model  of  Fay  and  Riddell  in  Refer¬ 
ence  5.  The  results  obtained  indicate  no  important  effects  for  catalytic 
or  near  catalytic  surfaces;  and  thus,  for  the  analysis  described  herein 
the  gas  is  treated  as  thermally  and  calorically  perfect. 

Stetson  in  Reference  6  has  presented  an  empirical  correlation  of  the 
results  of  the  numerical  solution  for  momentum  thickness,  0,  of  Kemp, 
Rose,  and  Detra  which  is  given  by  Equation  24  (Appendix  I).  This  equa¬ 
tion  applies  to  two-dimensional  or  axisymrnetric  flows  and  can  be  used 
with  reasonable  accuracy  for  quasi -symmetric  flows  such  as  that  over 
elliptical  cones.  Knowing  the  equation  for  the  momentum  thickness,  0, 
and  knowing  the  inviscid  flow  properties  at  the  edge  of  the  boundary 
layer  (from  Item  3  below),  tho  momentum  thickness  and  the  momentum 
defect 


M  =  f  p  u2Bd£ 

I  g  0 

''circumference 

can  be  evaluated.  The  details  of  this  development  are  given  in  Appendix  1. 
b.  Turbulent  Boundary  Layer 

The  state  of  knowledge  of  turbulent  boundary  layers  is  not  as  well  deve¬ 
loped  as  that  of  the  laminar  case.  Although  there  have  been  many  ap¬ 
proximate  treatments  of  turbulent  flow,  most  of  which  can  be  supported 
to  one  degree  or  another  by  some  of  the  existing  experimental  data, 
there  are  no  "exact"  theoretical  solutions  such  as  that  of  Kemp,  Rose, 
and  Detra  for  laminar  flow. 

One  approximate  solution,  the  results  of  which  lend  themselves  to  the 
present  problem,  is  that  of  Reshotko  and  Tucker  (Reference  7).  This 
solution  uses  the  momentum  integral  and  moment -of-momentum  equations 
and  is  simplified  by  using  Stewartson  s  transformation  from  Reference  8. 
In  order  to  solve  these  two  equations,  a  skin  friction  relation  must  be 
used.  Reshotko  and  Tucker  chose  the  Ludwieg -Tillman  relation  (Refer¬ 
ence  9)  in  a  form  suitable  for  compressible  flow  with  heat  transfer 
through  application  of  the  reference  enthalpy  concept.  The  above  equations 
were  further  simplified  through  the  use  of  an  approximate  shear  stress 
distribution  and  the  power  law  velocity  profile. 

The  moment -of-momentum  equation  is  needed  in  order  to  account  for 
pressure  gradient  effects  on  the  boundary  layer  velocity  profile.  The 
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method  o;  Reshotko  and  Tucker,  a9  applied  to  insulated  surfaces,  is 
quite  well  founded.  However,  for  noninsulated  or  nonadiabatic  wall 
cases,  the  method  (though  qualitatively  correct)  is  hased  on  some  specu¬ 
lative  assumpti'  ns.  It  is  still  anticipated,  however,  that  for  such  cases 
the  method  will  yield  reasonable  quantitative  results.  The  method  cer¬ 
tainly  is  representative  of  the  best  solution  that  can  be  obtained  within 
the  present  status  of  turbulent  flow  theory. 

Reshotko  am’  Tucker  have  presented  their  results  in  terms  of  a  trans¬ 
formed  momentum  thickness,  0tr,  given  by  the  equation 

e  1  V  “  H-1  M«2'3  • 

If  the  boundary  layer  is  turbulent  over  the  entire  body,  the  solution  be¬ 
gins  at  the  stagnation  point.  However,  if  the  boundary  layer  is  initially 
laminar,  the  calculation  must  begin  at  the  transition  point.  The  transi¬ 
tion  point  can  be  specified  either  in  terms  of  a  location  or  a  transition 
Reynolds  number,  Reg,  or  it  can  be  determined  from  the  transition  cri¬ 
terion  developed  u«  Item  3,  b.  Thi3  criterion  has  been  prog  i  aiuinsd 
along  with  the  turbulent  boundary  layer  solution.  The  details  of  the  solu¬ 
tion  for  boundary  layer  momentum  thickness  and  momentum  defect  are 
given  in  Appendix  I  along  with  the  basic  elements  of  the  programmed 
boundary  layer  solution. 

3.  SUPPORTING  ANALYSES 

a.  Forebody  Inviscid  Flow  Field 

In  order  to  evaluate  the  boundary  layer  momentum  thickness  and  momen¬ 
tum  defect,  the  inviscid  flow  field  characteristics  at  the  edge  of  the 
boundary  layer  must  be  determined.  In  particular,  for  the  programmed 
boundary  layer  analysis,  the  following  flow  properties  must  be  evaluated 
along  the  foreoody:  (1)  Pe/PD,  (2)  P0/Ps,  and  (3)  Me.  A  general  com¬ 
putation  procedure  was  developed  to  calculate  the  flow  properties  along 
the  forebody  based  only  on  inputs  of  free  stream  conditions  and  forebody 
geometry. 

The  procedure  is  applicable  for  a  wide  range  of  free  stream  conditions 
and  is  suitable  for  quasi -symmetric  geometries  as  weii  as  for  axisyrr. - 
metric  and  two-dimensional  bodies.  A  tangent  cone  method  of  solution 
was  used  to  obtain  the  flow  properties  for  sharp-nosed  bodies  and  a 
combination  of  Newtonian  and  tangent  cone  methods  was  used  for  a  blunt- 
nosed  configuration.  To  this  end,  a  Taylor -Maccoll  conical  flow  solution 
was  programmed  along  with  a  logic  procedure  for  the  application  of  the 
tangent  cone  and  Newtonian  methods  to  quasi -symmetric  forebody  geo¬ 
metries.  The  logic  procedure  is  shown  in  Appendix  II. 

The  use  of  Newtonian  theory  and  tangent  cone  theory  is  sufficient  to 
yield  solutions  for  pressure  and  Mach  number  within  10  to  20  percent 
accuracy,  which  is  quite  sufficient  for  engineering  purposes.  Since  the 
momentum  defect  varies  as  the  square  root  of  the  pressure,  a  20-per¬ 
cent  error  in  pressure  would  cause  less  than  a  10-percent  error  in 
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momentum  defect.  Note  also  that  the  effects  of  boundary  layer  separa¬ 
tion  on  the  pressure  are  not  considered  in  this  analysis. 

b.  Boundary  Laver  Transition 

In  order  to  make  use  of  the  turbulent  boundary  layer  analysis,  it  is 
necessary  to  determine  the  location  at  which  transition  occurs.  As  a 
part  of  the  computer  program  developed  to  predict  boundary  layer  mo¬ 
mentum  defect,  three  choices  exist  for  determining  transition: 

1.  Input  the  transition  point  location 

2.  Input  the  transition  Reynolds  number,  Regtr 

3.  Compute  R-e0tr  based  on  a  programmed  transition  criterion 

The  criterion  developed  to  predict  transition  was  based  on  a  correlation 
by  Dr.  R.  Nerem,  consultant,  of  the  form 

Rep  «  1000(1  -  0.  12  M2  +  0.023  M3)1/2 
tr 

(-2.29  +  17.  38  FT  -  18  K2  +  3.91  H3)1^2  ,  (1) 


where 


h  =  h  /h 
w 


aw 


For  regions  of  adverse  pressure  gradients,  such  as  in  the  vicinity  of  a 
compression  surface,  the  Reg  is  replaced  with  the  critical  Reynolds 
number  given  by  r 


Reg  =  163  (1  -  0.  12  M2  +  0.023  M3)1//2 

critical 

(-2.29  +  17.  38  IT  -  18  H2  +  3.91  KJ)  1/2  .  \2) 

These  equations  have  been  programmed  as  part  o'.  the  boundary  layer 
calculation  procedure  and  are  used  if  transition  location  or  transition 
Re0  is  not  specified.  The  evolution  of  the  above  equations  and  a  general 
discussion  of  boundary  layer  transition  may  be  found  in  Appendix  III. 
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SECTION  III 
WAKE  ANALYSIS 


1.  WAKE  REGIONS 

Analysis  of  thi  characteristics  of  high-speed  *  akes  begins  by  examining 
the  various  regions  that  comprise  the  conventional  high-speed  wake. 

The  five  dirtinct  regions  of  the  high-speed  w  ke  shown  in  Figure  1  are: 

1.  Recirculating  base  flow 

2.  Free  shear  layer 

3.  Rear  stagnation  point  and  the  neck  region 

4.  Inner  viscous  wake  (downstream  of  the  neck) 

5.  Outer  inviacid  wake 

Each  of  these  regions  has  its  own  peculiar  characteristics.  The  recir¬ 
culating  base  flow  region,  for  example,  ii»  boundary  layer  fluid  entrained 
due  to  the  adverse  pressure  gradient  associated  with  the  deflection  of 
fluid  at  the  axis  of  symmetry  back  into  the  streamwise  direction.  The 
free  shear  layer  is  also  boundary  layer  fluid;  this  fluid  initiates  the  de¬ 
velopment  of  the  neck  region,  which  is  characterized  by  steep  gradients 
due  to  the  presence  of  the  T'  empress  ion  shock. 

Aft  of  the  neck  region,  the  viscous  wake  develops  in  a  rather  orderly  and 
well-behaved  manner.  This  viscous  wake  is  the  subject  of  Item  3  of  thi«i 
section  in  which  the  conventional  boundary  layer  equations,  including  tl.e 
various  integral  relations,  are  applied  to  the  viscous  wake  and  a  momen¬ 
tum  integral  method  is  used  in  obtaining  a  solution 

The  outer  inviscid  wake  also  is  shown  in  Figure  1.  It  surrounds  the 
viscous  wake  and  reflects  the  nature  of  the  forebody  inviscid  flow  field 
and,  in  particular,  the  curvature  of  the  forebody  bow  shock  wave  and 
the  entropy  gradient  there  produced,  IN-*  inviscid  wake  i6  discussed  in 
detail  in  Item  4  of  this  section. 

2.  VISCOUS  WAKE  MODELS 

Considering  the  viscous  wake  in  more  detail  and  the  velocity  profiles 
shown  schematically  in  Figure  1  in  particular,  it  may  be  seen  that  for 
the  'onventional  high-speed  wake,  the  centerl.ne  X -direction  velocity, 

0,  increases  with  increasing  distance  lownstream  of  the  rear  stagnation 
point.  In  the  far  wake  [X/D  >>1),  the  centerline  velocity  actually  ap¬ 
proaches  asymptotically  to  the  lo  .1  free -stream  velocity.  For  this 
problem,  the  radial  velocity  gradient  is  zero  at  the  centerline  tor  all 
of  X,  as  are  also  the  gradients  in  temperature,  density,  and  the  otl.  -r 
fluid  properties.  On  the  other  hand,  the  outer  edge  of  this  viscous  wake, 
at  least  where  the  viscous  and  inviscid  wakes  are  clearly  separated 
from  one  another,  is  defined  in  the  usual  manner  for  boundary  layer 
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type  flows.  Here  the  properties  are  equal  to  their  inviscid  values,  and 
the  gradients  in  these  properties  are  zero.  Furthermore,  for  the  case 
where  viscous  and  inviscid  wakes  are  clearly  separated  from  one  another 
and  in  the  absence  of  any  pressure  gradient,  the  momentum  defect  will 
remain  constant  as  the  wake  develops  in  the  downstream  direction. 

When  a  decelerator  is  attached  to  a  forebody  by  a  cable  and  is  positioned 
in  such  a  conventional  high-speed  wake,  consideration  must  be  given  to 
the  extent  to  which  the  wake  characteristics  are  altered  due  to  the  inter¬ 
action  of  the  forebody  wake,  the  decelerator,  and  the  tow  cable  that  links 
these  two  bodies.  If  the  decelerator  was  located  in  the  immediate  neigh¬ 
borhood  o!  the  base  of  the  forebodv,  a  conventional  high-speed  wake 
would  not  have  the  opportunity  to  develop.  In  this  case  even  the  base 
flow  and  free -shear  layer  would  be  completely  altered  (see,  for  example, 
Figures  13  and  14  of  Reference  10),  and  the  proximity  of  the  decelerator 
to  the  base  of  the  forebody  might  result  in  the  presence  of  a  strong  ad¬ 
verse  pressure  gradient  in  the  largely  subsonic  cavity  region  separating 
the  two  bodies.  In  this  case,  the  wake  growth  is  explosive  and  the  de- 
ce'erator,  being  immersed  in  a  subsonic  cavity  type  of  flow,  produces  a 
low  drag  coefficient. 

If,  however,  the  decelerator  was  located  far  downstream  of  the  base  of 
the  forebody,  then  very  little  alteration  oi  the  wake  flow  would  be  ex¬ 
pected  di_’  to  the  interaction  of  the  decelerator  with  the  wake  itself.  In 
this  case  the  only  alteration  to  the  wake  (except  in  the  immediate  vicinity 
of  the  decelerator)  due  to  the  presence  of  the  tow  cable  or  riser  line. 

If  it  is  assumed  for  the  limiting  case  described  above  that  the  cab’e  dia¬ 
meter  is  very  small  compared  with  the  wake  diameter,  then  the  amount 
of  work  done  on  the  wake  due  to  skin  friction  on  the  cabie  is  neghgiDie. 
For  this  case,  in  the  absence  cf  any  pressure  gradient,  both  the  wake 
momentum  defect  and  momentum  thickness  are  constant  in  the  axial 
direction,  exactly  as  in  the  conventional  high-speed  wake.  However,  at 
least  for  laminar  flow  here,  the  similarity  be  twee'-,  these  wakes  ends 
because  m  the  present  case,  1 .  e ■  for  the  wake  with  a  trailing  decelerator 
in  it,  the  tow  cabie  forces  the  wake  centerline  velocity  to  be  zero.  As¬ 
suming  further  that  the  wake  profiles  are  similar  at  all  axia1  stations, 
then  the  ratio  of  the  wake  momentum  thickness  to  *he  wake  diameter  is 
the  same  at  all  axial  stations,  and  the  conclusion  is  that  the  wake  dia¬ 
meter  is  constant  with  axiai  position.  In  other  words,  for  this  particular 
case  there  is  no  wake  growth  such  as  is  present  in  the  conventional  high¬ 
speed  wake.  The  wake  is  cylmdrically  shaped  and  docs  not  change 
characteristics  in  the  axiai  direction.  This  is  illustrated  in  Figure  2. 

For  the  turbulent  viscous  wake,  however,  this  cylindrical  wake  model  is 
not  applicable.  The  turbulent  v  ake  in  the  presence  of  a  tow  cable  or 
riser  line  must  include  both  a  laminar  suoiayer  and  an  outer  turbulent 
mixing  layer  In  this  case,  the  efiects  of  the  tow  cable  are  solely  limited 
to  the  sublayer  region.  This  region,  though  not  present  in  the  conven¬ 
tional  wake,  is  such  that  it  allows  the  no-»hp  condition  to  be  satisfied 
at  the  wall.  i.  e .  the  cable  surface.  At  this  region’s  outer  edge,  the 
velocity  is  eq”al  to  that  of  the  turbulent  outer  layer.  This  outer  layer, 
on  she  other  hand,  not  being  constrained  by  the  nc-slip  condition  at  the 
centerline,  grows  as  in  the  case  for  the  conventional  turbulent  wake. 
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Figure  2  -  " Cylindrical "  Wake  Model 
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For  the  condition  of  a  contant  viscous  wake  momentum  defect  (the 
assumption  of  negligible  work  due  to  cable  skin  friction  is  still 
made  here),  the  velocity  at  the  inner  edge  of  the  turbulent  layer 
increases  just  as  the  centerline  velocity  increases  for  the  conven¬ 
tional  turbulent  wake. 

During  the  first  series  of  wake  survey  tests  (WT-IIA)  conducted  under 
the  EUREKA  program,  .•  ake  profile  measurements  were  obtained  behind 
an  axisyrmmetric  cone -cylinder -flare -cyl>  ,de r  forebody  with  and  without 
a  simulated  riser  line-  These  measurements  are  summarised  in  Fig¬ 
ures  3  and  4  where  impact  pressure  t  rofiles  at  four  different  X/D  sta¬ 
tions  are  presented  for  Mach  numbers  of  2  and  5.  There  is  virtually 
no  difference  between  the  properties  o  1  these  two  kinds  of  wakes.  Fur¬ 
thermore,  there  is  an  obvious  wake  growth  and  the  centerline  region 
veloc*  /  may  be  show,  to  increase  with  X/D  for  both  cases.  The  turbu¬ 
lent  wake  in  the  presence  of  a  riser  line  is  illustrated  in  Figure  5. 

The  discussions  presented  here  concerning  the  nat  re  of  both  laminar 
and  turbulent  wakes  in  the  presence  of  a  tow  cable  or  riser  line  are  based 
on  only  the  preceding  arguments  and  an  extremely  limited  amount  of  ex¬ 
perimental  data-  In  spite  of  this,  these  models  should  prove  useful  n 
engineering  calculations.  They  do  appear  to  be  applicable  for  a  wide 
range  of  conditions  in  terms  of  Reynolds  number,  Mach  number,  etc. 

They  are,  'owever,  limited  to  the  case  where  the  radius  of  the  tow  cable 
or  riser  line  u  much  less  than  the  wake  momentum  thickness  and  to 
X/D  >  3,  avoiding  the  "strong  interaction  region." 

3.  VISCOUS  WAKE  ANALYSIS 

a.  Wake  Flow  Assumptions 

An  integral  method  of  solution  was  applied  to  the  problem  of  determining 
the  properties  of  both  laminar  and  turbulent  viscous  wakes  with  parti¬ 
cular  attention  to  the  problem  of  the  wake  flow  preceding  a  trailing  aero¬ 
dynamic  decelerator.  In  the  paragraphs  that  follow,  the  theory  will  be 
summarized  for  the  different  types  of  wake  flows.  The  important  assum¬ 
ptions  contained  in  the  theory  are  as  follows: 

1.  The  fluid  is  thermally  and  caiorically  perfect,  and 
there  are  no  chemical  reactions  or  other  real  gas 
phenomena. 

2.  The  Prandtl  number  is  unity  (pr  =  1). 

3.  Prandtl's  concept  of  viscous  How  phenomena  is  valid 
for  the  high-speed  compressible  wake  such  that  gra¬ 
dients  in  the  streamwise  or  axta!  direction  are  much 
smaller  than  those  normal  to  the  wake  axis  and 
boundary  layer  type  equations  may  be  used. 

The  details  o:  the  base  flow  and  tree  shear  layer 
regions  may  be  largely  ignored  Thus,  the  reg.on 
of  validity  for  the  present  analysis  must  be  con¬ 
sidered  to  extend  from  somewhere  in  the  vicinity 
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o£  the  wake  neck  or  further  downstream  and  on 
down  into  the  far  wake. 

5.  Effects  due  to  the  existence  of  an  external  pressure 
gradient  are  negligible. 

Note  that  this  last  assumption  is  not  intended  to  restrict  the  solution  to 
flows  with  zero  pressure  gradient.  The  assumption  only  implies  that 
.ocal  acceleration  in  the  inner  wake  itself  is  much  more  important  than 
acceleration  due  to  an  external  pressure  gradient.  This  assumption 
allows  the  momentum  integral  equation  to  be  simply  written  as 


fix 

A 


Am 


<»  6j)  =  Its 


pu  (Uj  -  u)  rdr 


(3) 


constant  with  respect  to  X. 

Here  9  is  the  momentum  thickness  and  Equation  3  *s  simply  a  statement 
to  the  effect  that  the  wake  momentum  defect  is  constant  in  the  axial  di¬ 
rection.  Since  the  pressure  gradient  hac  been  assumed  to  be  negligible, 
then  there  is  no  external  work  being  done  on  the  wake,  and  the  momen¬ 
tum  defect  must  be  constant. 

As  a  check  on  the  validity  of  the  assumption  of  constant  wake  momentum 
defect,  experimental  wake  momentum  defect  values  were  calculated 
from  wake  survey  data  obtained  during  wind  tunnel  tests  conducted  under 
this  program.  In  looking  at  the  viscous  wake,  it  was  found  that  the  width 
of  the  wake  decreased  with  increasing  Mach  number.  For  this  reason, 
the  viscous  wake  momentum  defect  comparisons  were  conducted  at 
Mach  2.  The  comparisons  were  made  for  an  axisymmctric  spherical 
nose  cone-cylinder -fiare  forebody  at  a  Mach  number  of  2  and  a  Reynolds 
number  of  0.  369  X  10^/in.  The  wake  momentum  defects  were  computed 
at  several  wake  stations  using  a  data  reduction  program,  which  reduced 
wake  survey  pressure  and  temperature  data  to  local  velocities.  The 
data  reduction  program  used  a  linear  variation  between  experimental 
data  points  in  computing  momentum  defects.  The  edge  of  tha  viscous 
wake  was  determined  by  using  the  criterion  that  at  the  edge  of  the  viscous 
wake,  the  total  pressure  would  be  equal  to  the  total  pressure  immediately 
behind  a  normal  forebody  bow  shock.  The  data  reduction  program  those 
the  data  point  having  a  total  pressure  closest  to  this  pressure.  The 
spacing  of  the  experimental  data  points  was  equal  to  one -tenth  of  the 
forebody  base  diameter.  The  experimental  wake  momentum  defects 
also  were  compared  with  forebody  momentum  defect  calculated  using 
the  method  described  in  Section  II.  The  results  of  the  comparison  are 
shown  in  Table  I. 

The  tabulated  values  indicate  that  the  method  used  for  predicting  the  edge 
of  the  viscous  wake  is  not  accurate  enough  to  produce  accurate  experi¬ 
mental  wake  momentum  defect  results.  The  X/D  =  3  data  are  too  near 
the  wake  neck  to  be  a  useful  profile  for  this  study.  The  general  decrease 
in  the  momentum  defect  values  from  X/D  =  5  to  X/D  =  11  is  the  result  of 
the  inaccuracy  in  selection  of  the  nearest  experimental  data  point  as  the 
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TABLE  I  -  VISCOUS  WAKE  MOMENTUM  DEFECT  COMPARISON 


M 

00 

He  ^ 
to /in. 

■* 

X/D 

Momentum 
defect  in 

Y  plane 
(lb) 

Wake 

radius, 

Y/D 

Momentum 
defect  in 

Z  plane 
(lb) 

Wake 

radius, 

Z/D 

Fcrebody 

momentum 

defect 

(lb) 

2 

0.  369  X  10b 

3 

5.68 

0.  5 

2.  98 

0.  5 

1. 203 

2 

0.  369  X  10b 

5 

4.46 

0.  5 

1. 65 

0.  5 

1.203 

1 

2 

0.  369  X  106 

7 

4.  10 

0.  5 

1.02 

0.4 

1. 203 

2 

0.  369  X  10b 

9 

2.29 

n.  5 

0.964 

0.4 

1. 203 

2 

0.  369  X  10b 

11 

1.  73 

0.  5  1 

0.  786 

0.4 

1.203 

viscous  wake  edge.  Since  the  wake  is  growing,  the  value  of  the  momen¬ 
tum  defect  within  a  given  circle  will  decrease  with  increasing  X/D.  The 
experimental  data  also  may  contain  errors  induced  by  the  forebody 
mourting  strut  that  spanned  the  wind  tunnel  in  the  Z  plane.  Within  the 
limits  of  the  accuracy  of  the  experimental  data  and  the  computation 
method,  the  momentum  defect  is  assumed  to  be  constant.  In  an  attempt 
to  improve  the  method  for  computing  the  momentum  defect  in  the  wake, 
a  polynomial  curve  fit  technique  was  used  to  smooth  the  experimental 
data  and  therefore  improve  the  accuracy  of  the  calculation.  Using  the 
curve  fit  technique  i.i  the  wake  momentum  defect  calculations  did  not 
have  an  appreciable  effect  on  the  results. 

Assumption  5  listed  above  thrs  does  not  limit  the  solution  to  wake  flows 
with  zero  pressure  gradient.  However,  the  validity  of  the  solution  cer¬ 
tainly  does  decrease  with  increasing  pressure  gradient.  For  blunt 
bodies  with  strong  pressure  gradients,  the  solution  should  be  considered 
as  being  highly  approximate. 

By  ignoring  the  details  of  the  bas.*;  flow  and  shear  layer  regions  and 
assuming  that  no  external  work  is  done  on  the  wake  fluid  in  these  regions, 
then  the  constant  in  Equation  3  may  be  evaluated  by  equating  the  momen¬ 
tum  defect  of  the  boundary  layer  at  the  base  of  the  forebody  to  the  mo¬ 
mentum  defect  of  the  wake.  This  may  be  stated  mathematically  as 

(77  0bD)  =  (77  0^)  ,  (4) 

vhere  the  subscript  1  refers  to  the  properties  associated  wit’.,  the  outer 
edge  of  the  wake  and  the  subscript  B,  to  properties  at  the  base  of  the 
vehicle  where  the  boundary  layer  separates.  The  symbols  Pg  and  ug 
represent  inviscid  properties,  and  D  is  the  forebody  bast,  diameter. 

Next,  the  Dorodnitzyn  transformation  is  introduced  in  order  to  obtain  a 
solution  to  the  axial  direction  momentum  equation.  In  the  integral 
method  of  solution,  it  is  ..ssumed  that  the  velocity  profile  in  the  viscous 
wake  may  be  expressed  as  a  fourth-order  polynomial  of  the  form 
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Applying  the  appropriate  boundary  conditions  to  the  above  equation,  a 
solution  for  the  viscous  wake  velocity  profile  in  the  transformed  coordi¬ 
nate  system  may  be  obtained.  Considering  the  energy  equation  for  the 
assumption  of  Pr  =  1  and  applying  the  appropriate  boundary  conditions 
to  the  energy  equation,  the  enthalpy  profile  in  the  transformed  coordinate 
system  is  obtained  as  a  funct  on  of  £  the  form  parameter  as  derived  in 
Appendix  IV. 

The  details  of  the  above  development  are  given  in  Appendix  IV  and  apply 
to  both  laminar  and  turbulent  viscous  wake  flows.  The  necessary  equa¬ 
tions  for  determining  the  properties  of  a  high-speed  compressible  wake 
were  derived  using  an  integral  technique  and  the  Dorodnitzyn  transfor¬ 
mation.  In  order  to  describe  the  entire  wake  flow  field,  the  inviscid 
flow  field  must  be  determined  (see  Item  ,  below)  and  that  an  estimate 
be  made  of  the  vehicle  boundary  layer  momentum  thickness  (Section  II 
above). 


b.  Laminar  Axisymmetric  Wake 


Appendix  IV  describes  the  case  of  both  the  conventional  laminar  wake 
and  the  cylindrical  laminar  wake.  In  the  appendix,  the  relationship  be¬ 
tween  £  and  x  is  determined.  For  each  case,  it  is  necessary  to  specify 
a  relation  for  the  viscosity  coefficient  and,  lor  both  laminar  flow  cases, 
the  viscosity  is  assumed  to  be  linearly  dependent  with  respect  to  fluid 
temperature;  i.  e.  ,  /v/RT  -  constant.  For  the  cylindrical  wake  case, 
uq/ui  =  0,  which  is  equivalent  to  £=24  (Appendix  IV).  Thus,  the  pro¬ 
perties  of  the  cylindrical  wake  are  solved  by  setting  £  =  24,and  £will 
not  depend  on  X  for  this  case  since  there  will  be  complete  similarity 
between  radial  profile  at  each  X  direction  station. 

Having  solved  for  <f  as  ^  function  of  X  (Appendix  IV/  together  with  the 
equations  derived  in  that  appendix  for  the  velocity  and  enthalpy  profiles 
and  the  relationship  between  the  transformed  and  physical  radial  coor¬ 
dinates,  the  properties  of  the  laminar  viscous  wake  are  completely 
known. 

Using  this  theory,  as  shown  in  Figure  6,  reason,  ble  agreement  is  at¬ 
tainable  even  in  the  near  wake  region  where  the  j.  esent  theory  should 
give  the  least  reliable  results.  A  further  check  on  the  laminar  wake 
theory  is  presented  in  Figure  7  where  the  wake  growth  behind  a  conical 
body  is  presented  as  a  function  of  distance  downstream  of  the  body. 

Figure  7  also  shows  both  ballistic  range  experimental  data  from  Refer¬ 
ence  1 1  and  the  result  of  a  theoretical  calculation  using  the  laminar  wake 
method  developed  in  Appendix  IV.  As  may  be  seen,  excellent  agreement 
exists . 

Unfortunately,  good  data  does  not  exist  at  the  present  lime  for  the  laminar 
wake  in  the  presence  of  a  accelerator  with  which  the  laminar  cylindrical 
wake  results  may  be  compared.  Note  that  during  the  wind  tunnel  tests 
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Figure  7  -  Comparison  of  the  Momentum  Integral  Method  Solution  with 
Experimental  Laminar  Far  Wake  Data 


performed  as  a  part  of  this  program  only  turbulent  wakes  were  obtained. 
It  is  felt,  however,  that  since  the  conventional  laminar  wake  results  in 
Figure  6  are  for  conditions  similar  to  these  of  the  cylindrical  wake  (for 
example,  the  centerline  velocity  is  low  in  the  near  wake  region  and  the 
velocity  profile  is  similar  to  that  of  a  cylindrical  wake);  and  since 
reasonably  good  agreement  exists,  then  the  laminar  cylindrical  wake 
results  should  lead  to  equally  good  agreement. 

Turbulent  Axisymmetric  Wake 


For  the  turbulent  wake  case,  a  form  for  the  viscosity  relationship  must 
be  chosen  different  from  that  used  for  the  laminar  case.  In  turbulent 
flow,  the  viscosity  coefficient  is  only  an  effective  one  (artificially  intro¬ 
duced  due  to  the  manner  in  which  the  governing  equations  are  expressed) 
and  is  related  to  the  momentum  transfer  caused  by  the  fluctuating,  mixing 
type  motion.  Many  different  mathematical  forms  for  the  turbulent  visco¬ 
sity  coefficient  have  been  suggested  in  the  literature,  and  tour  of  these,  as 
applied  to  the  wake  centerline,  are  summarized  below  (these  four  models 
are  virtually  identical  to  those  discussed  in  Reference  13): 

"0  =  K  'Vo  ■ 

A0  =  K  A  ( PjUi  - 
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and 


UQ  =  K  (Uj  -  Uq)  , 

*0  =  K  Ap0(ul  *  U0)‘ 


While  only  the  above  models  were  considered  in  this  investigation,  the 
solution  is  not  restricted  to  these  viscosity  models.  In  all  of  the  above 
models,  the  value  of  the  constant  K  is  required.  The  first  model  relates 
the  turbulent  viscosity  coefficient  to  the  centerline  momentum  flux,  the 
second  model  form  involves  a  momentum  difference,  and  the  third  and 
fourth  models  involve  only  a  velocity  difference.  These  latter  models 
correspond  to  the  accepted  "classical"  form  for  incompressible  flow. 
However,  the  second  model  also  reduces  to  a  velocity  difference  for  the 
incompressible  case.  The  first  model  obviously  does  not  reduce  to  the 
classical  form  for  incompressible  flow.  In  fact,  the  first  model  did  not 
originate  from  classical  viscous  flow  work,  but  from  studies  of  jet¬ 
mixing  phenomena.  All  of  the  above  modes  may  be  expressed  using 
the  general  form  written  below: 


. .  .  ,  O.m 

"o  =  KA'iui  <7^ 


.n 


'0m  U  0 


n  -  (-!)“(-) 


(6) 


In  Table  II,  the  values  for  the  constants  m,  n,  and  p  are  presented.  In 
the  computer  program  for  the  viscous  turbulent  wake,  the  specification 
of  the  model  number  directs  the  program  to  use  the  proper  values  of  m, 
n,  and  p. 


The  value  of  K  must  be  determined  from  experiments.  In  low -speed 
flows,  it  ic  conventional  to  define  ,  an^  thus  K,  through  introduction 
of  the  Townsend  Reynolds  number,  Rj,  where 
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u  )r 
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(7) 


TABLE  II  -  CONSTANTS  FOR  TURBULENT  VISCOSITY  RELATIONSHIP 


M  del 

"0 

m 

n 

p 

: 

Remarks 

F.rst 

K  A  *quq 

0 

Momentum  flux 

Se  5nd 

K  A  ('iui  ‘  Vo1 

i 

Mome ntum  diffe  renc e 

1  hird 

K  A  />l  ( u  j  -  u q) 

o  ! 

i  0 

0 

Velocity  difference 

Fourth 

K  A  *0  (“1  '  "0) 

i  | 

uL 

0 

Velocity  difference 
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Here,  r*  is  the  so-called  half-width  at  half -depth,  i.  e-  ,  the  radius  at 
which  the  velocity  is  exactly  intermediate  to  the  axial  and  edge  values, 
and  Ry  is  the  constant  to  be  determined  from  experiment.  For  the 
secona,  third,  and  fourth  models  and  in  the  limit  of  either  incompressi¬ 
ble  flow  or  an  "incompressible  like"  flow,  i.  e.  the  far  wake  where 
PQ  -  Pj,  then  K  and  Ry  are  related  by 

K  =  .  (8) 

RT 

For  the  polynomial  velocity  profile  of  Equation  65  of  Appendix  D,  it 
follows  that 

* 

r  ^ 

s  n  =■  o.  4  , 

and  thus 


K  a 


(9) 


Typical  experimental  values  for  R^.  are  tabulated  in  Table  III  together 
with  the  value  of  K  calculated  from  Equation  9.  These  values  of  K  may 
be  compared  with  those  of  Zakkay  and  Fox  in  Reference  13  where  the 
first  model  is  similar  to  heir  Model  1,  the  second  model  to  their  Mo¬ 
del  2,  the  third  model  to  their  Model  2A,  and  the  fourth  model  to  their 
Model  4  Based  on  the  analysis  of  wind-tunnel  data,  Zakkay  and  Fox 
concluded  that  a  value  of  K  =  0.  01  for  the  second  model  and  K  =  0.  04  - 
0.  06  for  third  and  fourth  models  represented  a  best  fit  to  their  data. 
Zakkay  a..d  Fox  also  suggested  a  value  of  0.02  for  use  with  the  first 
model.  In  project  EUREKA,  the  best  correlation  with  wind-tunnel  wake 
test  data  has  been  obtained  for  values  of  K  on  the  order  of  0.005  to  0  01. 
Details  of  the  comparison  of  the  theory  with  experimental  data  will  be 
given  later  in  this  section. 

By  substituting  Equation  6  for  the  viscous  wake  into  Equation  63,  Appen¬ 
dix  IV,  for  ,  the  form  parameter,  the  relation  between  <f and  X  can  be 
determined. 


TABLE  HI  -  TOWNSEND  REYNOLDS  NUMBER  IN  WAKE 


Wake  type 

RT  ! 

1 

K 

Axisymmetric,  incompressible 

14  1  ! 

,  -u; 

Two  -dimensional,  incospressibie 

12.0  ! 

0.033 

Two  -dime  nsiona  1,  com  pres  s  ibie 

10  0  ’ 

|  '■  U4J 

Axisymmetric,  compressible 

13.7  | 

|  0.02b 
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If  the  present  method  of  solution  and  the  turbub  nr  viscosity  coefficient 
models  included  in  the  analysis  are  accepted,  then  the  resulting  equations 
for 4;  as  a  function  X,  together  with  the  equations  derived  in  Appendix  IV 
for  the  velocity  and  enthalpy  profiles  and  the  relationship  between  the 
transformed  and  physical  radial  coordinates,  are  sufficient  to  completely 
determine  the  properties  of  the  turbulent  viscous  w^ke.  Again,  the  pro¬ 
perties  of  the  forebody  boundary  layer  and  at  the  edge  of  the  viscous 
wake  are  assumed  to  be  known  from  the  fcrebody  computer  program. 

Using  the  above  method  of  analysis,  the  turbulent  viscous  wake  charac¬ 
teristics  were  determined  and  comparisons  were  made  with  the  velocity 
profiles  obtained  from  wind-tunnel  test  data  obtained  in  the  wake  of  the 
Arapaho  C  vehicle  with  nose  cone  (Appendix  V).  These  comparisons 
were  based  on  the  iterated  wake  solution  desciibed  in  Item  5  of  this  sec¬ 
tion.  In  Figures  8  through  2  J,  the  wake  velocivy  profile.,  are  shown 
plotted  versus  the  nondimensicnal  distance  (Y/D  or  Z/D)  measured  from 
the  centerline  of  the  viscous  wake.  Comparisons  ^re  shown  for  a  range 
of  Mach  numbers  from  2  to  6  for  two  Reynolds  number  conditions  end 
for  stations  behind  the  base  of  the  fcrebody  from  3  to  11  calibers,  which, 
is  the  region  of  interest  for  decelerator  application. 

The  theory  is  shown  a&  a  sMid  line  and  -epresents  the  total  wake  de¬ 
scription,  i.  e.  ,  both  viscous  and  inviscid  po-tions.  The  inviscid  wake 
comparisons  are  further  discussed  in  Item  4  of  this  section.  The  ex¬ 
perimental  velocity  data  points  shown  were  calculated  from  the  wake 
wind  tunnel  pressure  me  uurements.  obtained  as  ^  part  of  this  program. 
The  disparity  between  the  V/D -plane  date  and  the  Z/D-plane  data  is  at 
least  partially  due  to  the  presence  of  the  .orebody  mounting  strut  in  the 
Z/D  plane,  'however,  this  strut  also  may  influence  the  Y  /D  data,  and 
thus  data  for  both  planes  is  shown.  Note  also  that  at  X/D  =  3,  there 
may  be  added  inaccurac.es  in  the  data  due  to  the  influence  of  the  rake 
used  to  obtain  pressure  readings. 

The  best  correlation  with  the  experimental  EUREKA  wake  data  has  been 
obtained  using  the  third  viscosity  model,  =  K  &  p  \  (u,  -  ucl  with  a 
viscosity  coefficient  K  =  0  005  This  K  value,  while  below  th*  nor¬ 
mally  used,  was  justified  by  considering  that  the  5  iution  had  oeer.  ite¬ 
rated  as  explained  in  Item  5  of  this  section  thus  had  changed  the 
boundary  layer  approach  that  originally  defined  i  e  range  of  acceptable 
viscosity  coefficient  values  Using  K  o  005  tor  ihe  third  model  (shown 
in  Figures  8  through  15),  agreement  with  experiment  to  well  within  10 
percent  has  been  obtained  in  almost  all  cases  over  the  range  of  Mach 
numbers  and  wake  location  considered  in  this  program  and  for  both 
Reynolds  number  conditions  considered.  For  X/D  >  3.  the  correlation 
of  the  centerline  velocities  is  to  within  5  percent  in  all  but  one  case. 

Correlation  with  the  experimental  data  has  been  obtained  using  other 
viscosity  models,  the  results  using  the  fourth  model  are  sKown  in  Fig¬ 
ures  16  through  19  for  K  =  0.  010  and  in  Figures  20  through  23  for  K  = 
0.005  in  order  to  show  the  effect  of  changing  viscosity  model  and/or 
viscosity  coefficient.  Using  the  fourth  model,  the  results  show  that  the 
trend  of  the  predicted  centerline  velocity  is  affected  by  Mach  number. 

At  Mach  2  the  predicted  centerline  velocity  is  significantly  higher  than 
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Figure  11  -  Comparison  i ,i  Experimental  and  Calculaied  Wake  Profile 
(Mco^  5,  Pt  =  30.  5  psi,  Model  3,  K  =  0.005) 
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F’gure  16  -  Comparison  of  Experimental  and  Calculated  Wake  Profile 
(Mqq-  2,  Pt  =  2.8  psi,  Model  4,  K  =  0.01) 
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Figure  17  -  Comparison  of  Experimental  and  Calculated  Wake  Profiles 

(Mqo  =  3*  Pt  =  5.  9  psi,  Model  4,  K  =  0.01) 
oo 
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the  experimental  value.  As  the  Mach  number  is  increased,  the  degree 
of  correlation  changes  until  at  =  5  the  theory  significantly  under - 
predicts  the  velocity.  This  effect  is  not  seen  with  the  choice  of  the  third 
model.  For  the  same  K  value,  the  use  of  the  fourth  model  will  yield  a 
lower  centerline  velocity  than  will  the  third  model.  The  second  model 
for  which  no  data  are  shown,  will  result  in  even  higher  predicted  center¬ 
line  velocities  than  the  third  model.  The  first  model  was  excluded  due 
to  the  fact  that  it  deviated  from  the  proper  variation  of  velocity  with 
X/D  as  exhibited  by  the  other  models  (see  Appendix  IV).  With  each 
model  the  choice  of  the  viscosity  coefficient  K  will  influence  the  resultant 
velocity  profile.  As  K  is  increased,  the  velocity  predicted  by  a  given 
model  will  increase  and,  as  K  is  decreased,  the  predicted  velocity  will 
decrease. 

Prior  to  the  use  of  the  iterated  solution,  reasonable  correlations  between 
experiment  and  theory  were  obtained  for  a  value  of  K  =  0.  025  with  the 
third  and  fourth  models. 


d.  Wake  Transition 

In  Items  3,  ij  and  3,  c,  an  integral  method  for  calculating  the  properties 
of  both  laminar  and  turbulent  wakes  was  presented.  However,  up  to  this 
point  no  mention  has  been  made  of  the  conditions  under  which  the  viscous 
wake  would  be  laminar  or  turbulent.  Obviously,  if  the  forebedy  boundary 
layer  is  turbulent,  then  the  wcke  also  wiil  be  turbulent,  and  a  calculation 
of  the  wake  profiles  as  a  function  of  x/D  proceeds  as  described  in  Item 
3,  c.  However,  if  the  borebody  boundary  layer  is  laminar,  then  the 
viscous  wake  may  be  either  laminar  or  turbulent,  depending  on  whether 
or  not  the  flow  conditions  <*re  suitable  for  transition  to  occur. 

The  phenomena  of  transition  in  hypersonic  wakes  has  been  studied  ex¬ 
tensively  in  the  laboratory.  As  pointed  out  in  Reference  14,  the 
correlation  of  wake  transition  data  that  appears  to  indicate  the  most 
physical  insight  and  thus  be  the  most  general  is  that  of  Pallonc,  Erdcs, 
and  Eckerrnan.  In  their  work  they  used  a  transition  Reynolds  number 
based  on  local  properties  at  the  edge  of  the  wake,  a  velocity  equal  to 
that  at  the  edge  relative  to  the  centerline  velocity,  uj  -  uy,  and  on  the 
wake  diameter  at  transition,  The  correlati°n  °1  experimental 

data  that  they  present  is  thus  in  terms  of  the  transition  Reynolds  number, 
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and  the  effective  Mach  number, 
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Figure  24  shot  s  a  summarization  of  experimental  measurements  of  wake 
transition  Reynolds  number.  There  is  a  reasonable  correlation  of  the 
differert  experiments  and  the  wake  transition  may  be  predicted  to  within 
a  factor  of  two  in  Reynolds  number. 
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Figure  24  -  Experimental  Measurements  of  Wake 
Transition  Reynolds  Number 
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In  order  to  facilitate  calculations  using  the  correlation  of  Pallone,  Erdos, 
and  Eckerman,  an  approximate  curve  has  been  fitted  to  the  data  in  Fig¬ 
ure  24.  This  approximation  is 


Re 


=  1700  Mj' 


(12) 


TR 


This  curve  also  is  6hown  in  Figure  24,  and  reasonable  agreement  may 
be  see  to  exist. 

In  applying  Equation  12  to  a  laminar  cylindrical  wake,  uQ  will  be  zero 
due  to  no -slip  condition  on  the  tow  cable  or  riser  line.  Thus,  Re  p 
and  !&j  become  a  somewhat  more  conventional  "looking"  Reynolds  num¬ 
ber  ana  Mach  number  based  on  local  edge  conditions.  Furthermore, 
for  a  constant  pressure  wake  region,  Re  and  will  not  vary 

with  X/D.  The  viscous  wake  in  the  presence  of  a  riser  line  or  tow  cable, 
attaching  an  aerodynamic  decelerator,  will  be  either  wholly  laminar  or 
wholly  turbulent. 


Near  Wake  Characteristics 


The  general  characteristics  of  the  near  wake  were  illustrated  earlier 
in  Figure  1.  This  region  is  made  up  of  the  free  shear  layer,  the  base 
region  recirculating  flow,  and  the  beginnings  of  the  viscous  inner  wake 
at  the  rear  stagnation  point  and  on  down  through  the  neck;  when  con¬ 
sidered  in  total,  it  is  a  rather  complex  fluid  mechanic  phenomena. 


As  Keeves  and  Lees  have  pointed  out,  and  also  provided  in  Reference  15, 
a  relatively  clear  picture  of  the  structure  of  the  near  wake  now  exists 
as  a  result  of  the  many  experimental  studies  that  have  been  reported  in 
the  literature  (see,  for  example,  Reference  16  through  19).  However, 
theoretically,  the  state  of  the  art  has  not  progressed  as  far.  As  has 
been  pointed  out  by  Webb,  et  a,  in  Reference  20,  the  near  wake  is 
characterized  by  two  length  scales  Corresponding  to  the  existence  of 
separate  regions  of  high  shear  and  of  low  velocity.  Because  of  this  as¬ 
pect  of  the  problem,  the  resulting  governing  equations  are  nonsimilar 
in  fnrm  and  do  not  reduce  to  ordinary  differential  equations  when  pro¬ 
perly  transformed.  This  may  be  contrasted  wifh  the  far  wake  problem 
in  which  similar  solutions  are  obtainable. 


In  order  to  circumvent  this  complexity,  Lees  and  Reeves,  and  later 
Lees  working  with  Webb  and  others  at  TRW  Systems  (References  20 
through  22),  have  been  developing  multimoment  integral  methods  of 
solution  for  the  near  wake  problem.  Weiss  in  Reference  23,  on  the 
other  hand,  has  taken  an  approach  that  though  inherently  more  complex 
offers  the  hope  of  a  more  exact  solution  of  the  near  wake  problem. 

He  obtains  a  solution  by  dividing  the  near  wake  flow  into  three  regions: 
(1)  a  rotational  outer  flow  resulting  from  the  inviscid  expansion  of  the 
forebody  boundary  layer  at  the  base,  (2)  a  thin  viscous  "boundary"  layer 
above  the  dividing  streamline,  and  (3)  the  r circulation  region  itself. 

The  first  of  these  regions  is  solved  for  by  using  the  method  of  charac¬ 
teristics;  the  second,  using  a  modified  Oseen  solution  of  the  boundary 
layer  equations;  and  finally,  the  third  region,  by  a  finite  difference 


41 


AFFDL-TR -67-192 
Volume  II 


volution  to  the  full  Navi* r -Stokes  equations.  By  coupling  these  regions 
together  through  the  matching  of  boundary  conditions  at  their  interface* 
and  using  an  iterative  method,  complete  solutions  up  to  the  rear  stagna¬ 
tion  point  have  been  carried  out.  Although  this  method  has  not  yet  been 
extended  to  axisymmetric  bodies,  it  would  appear  to  be  the  type  of  ap¬ 
proach  that  could  lead  to  a  satisfactory  solution  of  the  near  wake  problem. 

For  the  present  investigation,  neglecting  the  n?ar  wake  details  by 
equating  the  momentum  defect  at  the  base  of  the  forebody  to  that  of  the  r 
viscous  wake  downstream  of  the  near  wake  region  appears  bo  be  a  satis*:  , 
factory  approximation  in  the  development  of  an  engineering  solution  for 
the  viscous  wake  characteristics  since  the  free  shear  layer  energy  trans¬ 
fer  to  the  recirculation  region  has  been  found  to  be  on  the  order  of  10 
percent  of  the  total  energy  transfer  from  the  bo  re  body  to  the  viscous 
wake.  However,  even  this  less  can  be  accounted  for  in  the  solution  by 
use  of  an  empirical  constant  as  will  be  shown  later  in  this  section. 

As  part  of  the  integral  wake  solution.  It  is  necessary  to  specify  the 
location  of  the  rear  stagnation  point  in  order  to  properly  relate  the  wake 
development  to  the  forebody. 

Although  the  rear  stagnation  point  has  been  chosen  in  this  instance  as 
the  starting  point  for  the  viscous  wake  solution,  any  other  suitable 
reference  location  <e.  g. )  wake  neck  or  sonic  point  may  be  used  if  suffi¬ 
cient  data  exist  to  specify  its  location,  Use  of  the  wake  neck  as  the 
starting  point  of  the  solution  could  yield  an  improvement  in  the  wake 
prediction  method  due  to  the  fact  that  by  starting  at  the  rear  stagnation 
point,  the  strong  gradients  in  the  wake  flow  are  ignored.  Further  ex¬ 
amination  of  the  starting  point  reference  should  be  considered  as  future 
program  improvements. 

In  Figure  25  most  of  the  available  data  in  the  literature  on  the  location 
of  the  rear  stagnation  point,  x0,  behind  bodies  at  supersonic  speeds  are 
shown.  Included  are  data  from  both  two-dimensional  and  axisymmetric 
bodies,  Mach  numbers  from  4  through  16,  and  for  a  wide  range  of  Rey¬ 
nolds  numbers.  Also  shown  is  an  approximate  fit  to  the  data  for  both 
the  laminar  and  turbulent  flow  regimes.  As  indicated  for  laminar  flow, 
the  rear  stagnation  point  moves  close;;  to  the  base  as  the  Reynolds  num¬ 
ber  increases.  This  movement  as  based  obviously  on  certain  allowance 
for  the  wide  data  scatter-  and  even  wider  discrepancy  for  a  few  individual 
points;  the  same  trend  also  is  noted  in  Reference  11. 

However,  this  movement  is  not  in  agreement  with  information  extracted 
from  available  theoretical  result*.  For  example.  Reeves  and  Lees 
predict  that  the  rear  stagnation  point  moves  away  from  the  body  with 
increasing  Reynolds  number  for  the  adiabatic  wall  case  (h^. *  Hj)  and 
the  cold  wail  case  (h#  s0,Z  H  j).  They  further  predict  that  the  rear 
•  {.agnation  point  will  be  closer  to  the  body  for  the  cold  wall  case  than  for 
the  adiabatic  wall  case  and  that  this  former  case  will  have  a  very  weak 
Reynolds  number  dependence-  This  picture  of  the  near  wake  perhaps 
can  be  reconciled  to  the  data  of  Figure  2§  by  noting  that  the  higher  data 
points  are  in  general  wind  tunnel  data  (e.  g. ,  Reference  17  and  22)  and 
tints  correspond  more  closely  to  the  adiabatic  wall  case,  while  ;  lower 
data  points  are  in  general  from  ballistic  range  and  shock -tunnel 
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experiments  and  closely  approximate  the  cold  wall  case.  This,  however, 
does  not  completely  explain  the  discrepancy  since  the  data  point  from 
Reference  18  also  was  obtained  in  a  wind  tunnel. 

Obviously,  the  "exact"  prediction  of  the  location  ci  the  rear  stagnation 
point  is  more  complicated  for  laminar  flow  than  presented  in  Figure  25. 
Still,  the  approximate  fit  in  Figure  25  should  prove  useful  for  engineer  - 
ing  calculations.  Furthermore,  for  cylindrical  wake  calculations,  the 
rear  stagnation  point  location,  xQ,  is  not  needed.  This  is  due  to  the 
complete  similarity  in  profiles  for  different  X/D  stations.  Also  Ehown 
in  Figure  25  are  measurements  of  the  location  of  the  rear  stagnation 
point  at  conditions  where  the  boundary  layer  at  the  base  of  the  forebody 
was  turbulent.  Although  there  are  only  a  limited  number  o*  measure¬ 
ments  for  this  case  of  turbulent,  flow,  an  approximate  fit. 


does  appear  to  correlate  well  with  the  data.  As  can  be  seen,  in  this  case 
there  is  no  apparent  dependence  on  Reynolds  number. 

As  indicated  earlier,  it  is  also  of  interest  to  compare  the  prediction  of 
the  rear  stagnation  point  enthalpy  with  available  experimental  data.  By 
combining  Equations  77  and  80  of  Appendix  IV,  the  expression  for  the 
rear  s  agnation  point  enthalpy  is 


£  i<  -  24}  =  1 


V 

H1  “1 


Equation  14  is  based  on  the  equating  of  the  momentum  and  total  enthalpy 
defects  of  the  wake  to  the  respective  quantities  of  the  forebody  boundary 
layer.  If  it  is  furthermore  assumed  that  ug  *  u j ,  the  simple  result  then 

is 


H  <£*  24)  h 
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(15) 


at  the  rear  stagnation  point.  :  Equation  15  is  compared  with  available  ex¬ 
perimental  data  for  laminar  flow  in  Figure  26.  The  measured  values  lie 
slightly  high*-  than  that  predicted  by  Equation  15. 

In  order  to  correct  the  present  theory  such  that  the  predicted  rear  stagna- 
trap  pfajpt  enthalpy  will  b-“  in  agreement  with  experiment.  Equation  4  and 
Equation  80/  Appendix  IV,  are  rewritten  as. 
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Figure  26  *  Expeiiment&l  Measurements  of  Rear  Stagnation 
Point  Enthalpy  for  Laminar  WaKe  Flow 


Here,  Kj  and  K7  are  Constanta  that  would  be  expected  to  be  greater  than 
or  equal  to  one  since  they  are  representative  -of  the  fractional  Increase 
in  (he  momentum  and  energy  defects  due  to  the  transport  Of  these  quanti¬ 
ties  aero**' .the  dividing  streamline.  Combining  the  above  two  equations 
with  Equation  7j,  it  follows  that 
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This  may  be  substituted  into  .-equation  77  with  the  result  that 
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where  K3  =  K^/KjUg. 

Equation  17  has  been  approximately  fitted  to  the  laminar  experimental 
data  using  a  value  of  K,  =0.9  and  is  shown  in  Figure  26.  Although  this 
value  of  Kj  is  Btrictly  based  on  experiment  and  has  no  real  theoretical 
basis,  its  use  in  Equation  17  for  laminar  wake  calculations  should  prove 
useful. 


For  turbulent  near  wake  uow,  there  is  only  a  smattering  of  data  that 
may  be  fourd  in  the  literature.  These  limited  number  01  data  points, 
however,  are  compared  with  Equations  i5  ?nd  17  for  both  =  0.9  and 
Kj  *  0.65,  in  Figure  27.  The  latter  value  of  does  give  agreement 
with  the  data-  The  limited  amount  of  data  prevents  any  real  recommenda¬ 
tion  for  the  use  of  =  0.65  for  turbulent  wake  calculations,  and  it  is 
obvious  that  more  data  and  more  comparisons  of  data  with  theory  are 
necessary  before  any  degree  of  confidence  can  be  placed  in  Equation  1? 
and  on  the  use  of  an  effective  K,  for  either  laminar  or  turbulent  flow. 

In  Appendix  IV  where  a  method  lor  viscous  wake  calculations  is  outlined, 
the  constant  is  left  as  an  input. 

As  was  noted  in  Item  3,  the  present  integral  theory  would  not  be  expected 
to  give  accurate  results  in  the  near  wake  region.  Thus,  the  lack  of 


Figure  27  -  Experimental  Measurements  oi  Rear  Stagnation 
Point  Enthalpy  for  Turbulent  Wake  Flow- 
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agreement  in  Figures  2f>  and  27  ie  not  surprising,  and  the  introduction 
of  an  empirical  constant,  K,,  into  Equation  17  only  represents  an  attempt 
to  improve  the  accuracy  of  me  present  «oluti  n  in  a  region  where  it  other¬ 
wise  woulu  not  be  expected  to  hold.  However,  in  spite  oi  the  failings  of 
the  p  csent  method  in  predicting  the  rear  stagnation  point  enthaipy,  torn* 
oariso.is  with  experimental  data  have  shown  that  in  the  prediction  of  neck 
diameter  for  a  laminar  wake  the  method  in  Item  3,  of  this  section  is 
quite  satisfactory. 

Before  leaving  this  discussion  of  the  near  wake,  one  last  comment 
should  be  made  with  regard  to  the  nature  of  the  flew.  It  has  been  found 
experimentally  that  the  near  wake  is  very  sensitive  to  the  base  geometry 
of  the  vehicle  and  also  to  the  injection  of  any  fluid  at  tin  base  into  the 
near  wake.  This  latter  effect  i6  obviously  of  concern  in  wind  tunnel  ex¬ 
periments  where  instrumentation  may  be  supporteQ  from  the  base  of  the 
forebody.  However,  it  also  opens  up  the  possibility  of  somewhat  control¬ 
ling  the  wake  characteristics  of  a  vehicle  through  fluid  injection.  It  is 
also  of  interest  in  the  decelerator  problem  where  due  to  the  manner  of 
decelerator  deployment  the  problem  of  interest  may  in  fact  be  one  in 
which  there  is  fluid  imection. 

In  Reference  24  T."..er  presents  measurements  that  show  a  50-percent 
change  in  the  location  of  the  rear  stagnation  point  due  to  a  nitrogen  in¬ 
jection  mass  flow  rate  that  is  only  2  percent  of  that  of  the  free  stream 
based  on  base  area  of  the  forebody.  Furthermore,  in  Referer.ee  25 
Bauer  shows  that  a  3-percent  injection  rate  using  helium  causes  the 
wake  recompression  shock  to  completely  disappear 

This  sensitivity  of  the  near  wake  to  almost  any  external  influence  casts 
doubt  on  the  applicability  of  any  wake  daU:  obtained  under  conditions 
where  componen’s  of  the  model  support  system,  e  g.,  struts,  wires,  or 
a  sting,  interfere  with  the  flow.  Thus,  future  efforts  should  include 
wind  tunnel  studies  where  such  interference  effects  are  not  present. 

4.  IN  VISCID  WAKE  ANALYSIS 

a.  Inviscid  Wake  Characteristics 

In  order  to  fully  describe  the  trailing  v,  ,»ke  behind  a  supersonic  body,  . 
both  the  viscous  inner  wake  and  -i:ivis/:id  outer  wake  characteristics, 
e.  g- ,  velocity,  density,  and  temperature  proxies,  must  t»c  known  The 
calculation  of  the  viscous  inner,  wake  depend  •»  ihe  forebody  boundary 
layer,  and  methods  lor  hand’u-ig  both  of  these  types  of  fluid  mechanic 
phenomena  have  been  previously  described  lor  both  -laminar  and  turbu¬ 
lent  flow  and  including  transition-  The  in  vise  ui  voter  wake  characteris¬ 
tics.  on  the  other  hand,  depend  on  fore  bod  v  bow  shock  wave  strength  a'ld 
shapi  and  the  subsequent  enpansion  over  thr  afi  portion  tv!  thr  body. 

b.  Analysis 


For  the  purpose  of  analysis,  the  wake  of  a  supersonic  vehicle  is  i  enr  i  - 
dered  where  it  :■  assumed  Uiat  the  vehicle  bow  shock  wave  shape  »nd 
the  pressure  distribution  along  the  wake  axis  **  known  from  -eithrr  ra¬ 
pe  rp  .ent  or  theory-  Fo'  blunt  bodscr-  at  hyperstjr.se  speed*,,  bias!  wave 
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theory  may  be  used  to  approximate  both  the  bow  shock  shape  and  the 
presaui  ,  distribution  on  the  wake  axis.  For  supersonic  speeds,  however, 
the  pressure  in  the  wake  may  oe  approximated  as  being  equal  to  P  ,  the 
tree  stream  pressure. 

From  conservation  of  mass,  a  mass  balance  may  be  written  that  relates 
the  mass  flow  in  the  free  stream  at  the  shock  through  a  disk  of  radius 
Rg  with  the  mass  flow  in  the  inviscid  wake  through  a  disk  region  of  radius 
r.  This  relation  i6 


V  ffR4 
oo  00  s 


purdr 


(18) 


The  integration  on  the  right-hand  side  of  Equation  18  starts  at  6  ,  the 
wake  displacement  thickness  instead  of  zero  in  order  to  account  for  the 
displacing  of  streamlines  by  the  viscous  wake  phenomena.  The  wake 
displacement  thickness  is  defined  by  the  relation 


f6 

*2  I 
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and  may  be  determined  using  the  viscous  wake  results. 

Equation  18  may  be  considered  to  be  a  relation  between  properties  on  a 
streamline  passing  through  a  poult  on  the  bow  shock  at  a  shock  radius 
Rg  and  properties  on  the  ume  streamline  downstream  of  the  body  at  a 
point  A,  a  distance  r  from  the  axis  of  the  wake.  With  the  bow  shock 
shape  known,  then  a  specific  value  of  Rs  corresponds  to  specific  proper  - 
ties  immediately  behind  the  bow  shock  for  the  streamline  passing  through 
that  point.  If  it  is  assumed  that  the  flow  is  isentropic  in  expanding  from 
the  shock  wave  on  downstream  to  point  A,  then  Equation  18  i-.  a  relation 
between  two  points  with  the  same  entropy,  or,  in  other  words,  the  same 
total  pressure.  The  assumption  of  isentropic  flow  behind  the  bow  shock 
means  that  any  secondary  shocks  (i.  e. ,  wake  recompression  shock  and 
iorebody  flare  shock)  produce  negligible  changes  in  entropy.  Although 
this  assumption  is  possibly  open  to  question,  it  simplifies  the  calculation 
of  the  inviscid  wake  properties  and  is  in  line  with  the  other  assumptions 
employed. 

Equation  18  thus  provides  a  means  for  calculating  the  total  pressure  or 
entropy  at  point  A,  a  distance  r  from  the  wake  axis.  Furthermore,  the 
total  enthalpy  is  known  since  the  inviscid  flow  may  be  considered  adiabatic. 
Fi  nally,  with  the  pressure  distribution  .Jong  the  axis  of  the  wake  known, 
an  estimate  of  the  static  pressure  in  the  inviscid  wake  can  be  made  by 
assuming  that  there  are  no  pressure  gradients  radially  in  the  wake  region 
i .  e ■  ,  dp/dr  -  0.  The  complete  state  of  the  gas  r.t  a  point  A  in  the  invis¬ 
cid  wake  thus  is  completely  speciiied  in  terms  of  the  total  pressure,  the 
total  enthalpy,  and  the  static  pressure  In  order  to  facilitate  calculations 
using  the  above  analysis,  the  following  transformation  is  introduced: 
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urdr  = 
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Equation  18  then  may  be  rewritten  as 
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where  ^  is  the  transxorrv  _^.  -placement  thickness,  and 
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Thus,  from  Equation  a  21  and  2?  and  assuming  <5  known  from  an  analysis 
of  the  viscous  wake,  the  entropy  and  total  pressure  at  a  point  A,  lying  a 
distance  n  fnrti  the  axis  i  i  *he  transformed  plane,  may  be  determined  as 
b^ing  equal  to  tha  immediately  behind  the  bow  shock  wave  at  a  radius 
R  .  'V.'.th  the  total  press  ire  or  entropy  known  for  the  streamline,  and 

wit',  an  e  .f’-nate  jf  he  static  pressure,  the  local  Mach  number  at  point 
A  may  h>  deterrr.in  d.  This,  together  with  the  total  enthalpy,  gives  the 
static  e  tiialpy  or  temperature.  The  complete  state  of  the  gas,  includ¬ 
ing  the  velocity  sol  density,  at  such  a  point  is  thus  specified,  and  the 
variation  in  fluid  properties  in  the  transformed  (n,  x^)  plane  may  thus 
be  determined.  In  order  to  transform  the  results  back  into  the  (X,  r) 
plane,  Equation  ?.0  may  be  used  with  the  result  that  r  and  n  are  related 
by  the  equation, 


rdr  = 


C  P  v 

I  00  CO 

X*  />’-* 


Since  Po^cd/p''1  known  as  a  function  of  n,  this  calculation  may  be 
readily  carried  out  using  numerical  techniques  for  the  integration.  The 
variation  of  such  fluid  properties  as  velocity,  Mach  number,  and  tem¬ 
perature  thus  may  be  calculated  as  a  function  of  r,  the  radial  distance 
from  the  axis  of  symmetry. 

Using  this  inviscid  wake  prediction  method  with  the  shock  shape  as  de¬ 
termined  from  wind  tunnel  tests  of  an  Arapahc  C  with  nose  cone  forebody 
as  an  input,  a  good  correlation  betv'een  theory  and  experimental  results 
has  been  obtained.  Figures  8  through  23  show  both  the  experimental 
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velocity  data  point*  obtained  by  reducing  the  measured  pressure  data 
and  the  inviscid  wake  theoretical  prediction*  for  Mach  numbers  of  ?.,  3, 

4,  and  5  for  two  Reynolds  number  conditions  and  for  X/D  values  between 
3  and  11.  The  starting  point  for  the  inviscid  wake  prediction  is  noted 
in  the  figures,  and  in  all  instances  the  theor.  is  within  10  percent  of  the 
experimental  values.  Thus,  the  degree  of  correlation  is  certainly  ac¬ 
ceptable  for  engineering  calculations. 

The  accuracy  of  the  inviscid  wake  theory  will  be  limited  to  some  extent 
by  the  accuracy  of  the  inputted  shock  shape.  In  performing  the  mass  flow 
balance,  theory  does  not  consider  the  embedded  shock  produced  by  the 
forebody  flare  or  the  recompressicn  shockwave  formed  in  the  wake  neck 
region  due  to  the  added  complexity  required  to  include  those  considera¬ 
tions  . 

Although  there  is  some  data  scatter  and  also  some  difference  between 
Y/D  plane  experimental  values  and  the  Z/D  plane  experimental  values 
due  in  part  to  the  fact  that  the  Z/D  values  were  measured  in  the  plane 
of  the  strut,  the  extent  of  the  correlation  between  theory  and  experiment 
is  sufficient  for  engineering  calculations  in  the  wake.  The  accuracy  of 
the  correlation  improves  for  larger  values  of  X/D.  In  the  region  of 
X/D  =  3  that  is  quite  near  the  wake  neck  and  the  recomprcssion  wake 
shock,  the  wake  prediction  method,  while  still  adequate,  is  Lss  accurate. 

5.  ITERATIVE  WAKE  SOLUTION 

In  order  to  completely  describe  the  characteristics  of  the  viscous  and 
inviscid  wake,  an  iterative  solution  was  performed  to  compute  the  viscous 
wake  profiles.  This  was  done  in  order  to  obtain  continuous  profiles  of 
the  wake  properties  acrocr  the  total  wa«  (viscous  +  inviscid)  since  the 
viscous  wake  edge  conditions  were  originally  calculated  based  on  the  flow 
through  the  normal  portion  of  the  bow  shock  wave  (as  per  boundary  layer 
analysis)  and  did  not  match  the  edge  value  obtained  frem  the  mass  flow 
balance  performed  in  the  inviscid  wake. 

The  iteration  procedure,  whicn  has  been  programmed,  computes  the 
viscous  wake  properties  basti  o*»  a  re /ised  value  ox  Mach  number,  Mj, 
at  the  edge  of  the  viscous  wake.  The  revised  value  of  M,  i3  obtained 
from  the  mass  flow  balance  performed  in  the  inviscid  wake  analysis  and 
is  that  value  of  Mi  predicted  by  the  mass  flow  balance  corresponding  to 
the  edge  location  predicted  by  the  initial  viscous  wake  run.  This  new  Mj 
value  is  used  as  an  input  for  a  second  pasa  (first  iteration)  through  the 
visco-  "  wake  program  and  the  inviscid  wake  program.  Although  the 
above  procedure  could  be  repea.eu  a*»  ii.ueliniie  number  oi  times,  initial 
results  have  shown  that  no  more  than  three  iterations  are  necessary  to 
obtain  complete  matching  of  the  viscous  wake -inviscid  wake,  edge. 

The  results  based  on  this  iterative  solution  are  shown  in  Figure*’  8 
through  23.  The  beat  correlation  for  these  turbulent  viscous  w:,.e  cases 
were  obtained  using  Model  3  with  a  viscosity  coefficient  oquai  to  0  005. 

A  correlation  for  laminar  viscous  wake  nas  not  been  made  using  the 
iterative  solivdon.  This  value  of  K  -  0.00a  is  below  the  ’•ang*  o'  values 
(0.02  to  0.06)  indicated  in  the  literature  as  shown  in  lt«.m  3,  c.  However, 
those  values  were  obtained  from  standard  boundary  layer  type  sol  tions. 
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By  performing  the  iteration  solution  that  change#  the  consideration  of  the 
flow  through  the  normal  shock  as  in  standard  boundary  layer  theory  to 
flow  through  a  locally  oblique  shock  as  used  on  this  program,  the  value 
of  K  yielding  a  satisfactory  correlation  has  been  changed.  Prior  to  the 
use  of  the  iterative  solution  with  the  viscous  wake,  reasonable  correla¬ 
tions  between  experiment  and  theory  were  obtained  for  value*  of  K  =  0.  025 
with  Models  3  and/or  4.  Until  mort  experimental  data  become  available 
on  the  effects  of  other  parameters  (such  <18  the  wake  neck  as  a  starting 
point  rather  than  rear  stagnation  point  mentioned  in  Item  3,  e_of  this 
section,  on  the  results  of  the  predictive  method,  an  acceptable  range  of 
K  values  cannot  be  exactly  determined.  Thus,  it  is  believed  that  the 
current  use  of  K  =  0.005  with  the  iterated  solution  is  acceptable  as  an 
engineering  solution  insofar  as  achieving  and  maintaining  a  general 
theoretical  approach  to  the  wake  description. 

FLOW  FIELD  COMPUTER  PROGRAM 

An  integrated  computer  program  has  been  developed  ami  programmed  in 
FORTRAN  IV  for  use  on  the  IBM  360-40  and  IBM  7094  computers.  The 
program  is  based  on  the  forebody  program  of  Section  II  and  the  wake 
analyses  of  Section  III. 

The  computer  program  can  determine  the  properties  of  the  inviscid  and 
viscous  wake  preceding  a  decelerator  and  consists  of  the  following  sub¬ 
programs: 

1.  Forebody  flow  field  (M^>  and  PjjJ 

2.  Boundary  layer,  laminar  and  turbulent  (0,  M*) 

3.  Viscous  wake,  laminar  and  turbulent  (M.  ,  P.  ,  and 

tl) 

4.  Inviscid  wake  (M^,  P^,  and  T^) 

5.  Iteration  (on  Items  3  and  4  above) 

The  following  inputs  must  be  specified  to  utilize  the  integrated  program: 

1.  Free  stream  conditions 

2.  Forebody  geometry  and  wall  temperature 

3.  Bow  shock  sV>ape 

4.  Boundary  layer  constant,  K  (0  or  1) 

5.  Wake  ccxilam.*,  K  and 


6.  Number  of  iterations  to  be  performed 

Optional  inputs  (wh'ch  will  be  calculated  in  the  program  if  not  specified) 
include: 
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1.  Boundary  layer  transition  Reynolds  number  or  transi¬ 
tion  point  location 

2.  Forebody  flow  field  (M^,  P^,  and 

Copies  of  this  computer  program  may  be  obtained  from  AFFDL(FDFR), 
WPAFB,  Ohio  45433  by  request. 
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SECTION  IV 

PARASONIC  PARACHUTE  PERFORMANCE 


1.  PERFORMANCE  FACTORS 

The  performance  of  PARASONIC  parachutes  is  described  according  to 
two  basic  criteria:  drag  and  stability.  The  primary  facers  that  influ¬ 
ence  the  performance  are  listed  below. 

1 .  Canopy  shape 

2.  Canopy  porosity 

3.  Free  stream  Mach  number 

4.  Free  stream  Reynolds  number 

5.  Forebody  shape 

6.  Downstream  location  of  the  canopy  inlet  (X/D) 

7.  Ratio  of  canopy  area  to  forebody  base  area  (A  /A^) 

The  first  two  parameters,  canopy  shape  and  orosity,  were  held  constant 
in  this  performance  evaluation.  All  parachute  testing  was  conducted 
with  PARASONIC  parachutes,  thereby  constraining  the  canopy  shape. 

The  selection  of  PARASONIC  parachutes  was  made  based  on  the  high 
degree  of  inflation  stability  and  sir"ctural  integrity  exhibited  by  these 
parachutes.  The  PARASONIC  parachutes  are  members  of  the  Hyperflo 
family  of  supersonic  parachutes  and  are  distinguished  from  standard 
Hype  rflo  parachutes  by  comparison  of  canopy  profiles.  The  PARASONIC 
parachutes  are  constructed  with  an  uninflated  shape  which  approximates 
the  fully  inflated  shape  of  a  standard  Hyperflo  parachute  (see  Figure  28). 
When  inflated,  the  PARASONIC  canopy  is  free  from  excess  material, 
thereby  reducing  "flagging"  loads  which  have  been  found  to  lead  to  pre¬ 
mature  canopy  failure.  The  porosity  was  held  constant  for  all  parachutes 
to  further  limit  the  number  of  variables. 

An  experimental  lust  program  was  conducted  to  determine  the  effects 
of  the  remaining  parameters  on  the  PARASONIC  parachute  performance. 
The  analysis  of  the  effects  of  these  parameters  is  presented  herein.  The 
experimental  test  program  is  discussed  in  detail  in  Appendix  V- 

2.  WIND-TUNNEL  TEST  SUMMARY 

The  PARASONiC  parachutes  described  above  were  deployed  behind  '  i<* 
three  forebody  configurations  shown  in  figure  29.  The  basic  Ar*t  .  io  C 
forebody  is  a  0.  l82 -scale  model  of  the  Arapaho  C  free -flight  test  vehicle 
and  consists  ol  a  probe  nose;  a  cylindrical  centerbody,  and  a  symmetrical, 
flared  afterbody.  The  Arapaho  C  with  nose  cone  forebody  has  a  center  - 
body  and  afterbody  identical  to  the  basic  Arapaho  C;  however,  the  probe 
nose  was  replaced  with  a  conical  nose  with  a  spherical  tip.  The  third 
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Figure  28  -  PARASONIC  Decelerator  Terminology 
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forebody  consisted  of  a  blunted  elliptical  cone  with  a  nose  shape  identical 
to  the  nose  used  on  the  second  forebody. 

Two  types  of  PARASONIC  canopies  were  tested.  One  type  was  con- 
structed  with  a  prewoven  Nome*  fine -mesh  cloth  as  the  roof  material, 
and  the  second  type  was  constructed  with  a  hand -woven  coarse -mesh 
roof.  Two  canopy  sixes  were  tested  so  that  a  canopy  area  to  forebody 
base  area  ratio  of  3.  425  could  be  obtained  for  ail  forebodies.  As  stated 
earlier,  all  canopies  were  constructed  with  the  same  total  aerodynamic 
porosity.  The  construction  details  for  the  parachute  models  are  shown 
in  Appendix  V. 

Two  series  of  parachute  performance  tests  we?e  conducted  under  this 
program.  The  first  series  or  tests  (WT-I1B)  was  of  sufficient  scope  to 
produce  enough  data  to  determine  the  effects  of  the  aforementioned  para¬ 
meters  on  the  parachute  performance.  The  second  series  of  performance 
tests  (WT-VI)  was  conducted  using  a  redesigned  drag  measurement  sys¬ 
tem  to  obtain  data  having  a  higher  degree  of  accuracy  as  the  comparative 
basis  for  the  drag  predictive  method  discussed  in  Section  V.  The  re¬ 
sults  of  this  series  of  tests  were  also  used  to  validate  the  results  of  the 
first  series  of  tests.  A  summary  of  the  test  conditions  for  both  series 
of  tests  is  shown  in  Table  IV. 

Results  obtained  from  the  parachute  performance  tests  are  presented  in 
Appendix  V.  Both  drag  and  stability  parameters  were  determined.  The 
base  pressure  to  free  stream  static  pressure  (Pb/P^)  values  were  re¬ 
corded  to  indicate  the  degree  of  modification  of  the  forebody  wake  due 
to  the  proximity  of  the  attached  decelerator.  Table  V  lists  the  P. /P 
ratios  for  th  three  forebodies  without  a  decelerator  immersed  in  the° 
inner  wake. 

Conpar^ion  of  the  values  in  Tables  V  with  the  results  presented  in  Ap¬ 
pendix  V  gives  an  indication  of  the  degree  of  wake  modification.  The  re¬ 
maining  parameters  listed  in  Appendix  V describe  the  stability  and  drag 
characteristics  of  the  parachute.  General  stability  comments  are  dis¬ 
cussed  ir.  lienrr  8,  below. 

3.'  MACH  NUMBER  EFFECTS 

Figures  30  and  31  illustrate  the  effects  of  the  tree  stream  Mach  number 
on  the  drag  coefficient  of  PARASONIC  models  behind  the  Arapaho  C  with 
nose  cone  and  the  blunted  elliptical  one  forebodies.  The  trend  of  de¬ 
creasing  drag  coefficient  with  increasing  Mach  number  shown  in  these 
figures  ft  typical  of  supersonic  parachute  performance.  The  Mach  num¬ 
ber  effset,  however,  is  interrelated  with  forebody  shapes,  parachute 
axial  location  (X/Dt.  and  canopy  to  forebody  base  area  ratio  (Ac/Afah 
The  major  interrelationship  of  interest  is  the  combination  of  Mach  num¬ 
ber,  X/D,  and  A  /A^  va’ues  at  which  lever,  modification  of  the  forebody 
wake  is  encountered,  li  nas  been  found  that  for  a  given  forebody  -para  - 
chute  combination,  the  axial  location  at  which  severe  modification  is 
encountered  increases  as  the  free  stream  Mach  number  increase*. 
Further  discussion  '  iteT  elation  of  the  free  stream  Mach  number 
and  other  influence  .eters  is  presented  in  applicable  sections. 
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TABLE  IV  -  PARACHUTE  PERFORMANCE  TEST 
CONDITIONS  SUMMARY 


Forebody 

Mesh 

— 

V\ 

M 

oo 

(pei) 

x/ 

Min 

Max 

Min 

WT-IIB 

Basic 

Small 

3.  425 

2,  3,  4,  5 

0.  5 

1.5 

4  5 

Arapaho  C 

Large 

3.  425 

2.  J.  4,  5 

0.  5 

1.5 

4.  5 

Small 

5.  9 

3.  4 

0.  5 

1.5 

6.0 

Large 

5.  9 

3.  4 

0.  5 

1.5 

6.  0 

Blunted  elliptical 

Small 

1. 99 

2,  3,  4.  5 

0.  5 

1.  5 

4.  0 

cone 

Large 

1. 99 

3,  4 

0  5 

1.  - 

4.0 

Small 

3.  425 

2.  3,  4.  5 

0.  5 

1.5 

4.  5 

Large 

3.  425 

2,  3.  4,  5 

0.5 

1.5 

4.  5 

Arapaho  C  with 

Small 

3.  425 

2.  3,  4.  5 

0.  5 

1.0 

4.  5 

nose  cone 

Large 

3.425 

2.  3,  4.  5 

0.  5 

1.0 

4.  5 

Small 

5.  9 

3.  4 

0.  5 

1.0 

6.  0 

Large 

5.  9 

3,  4 

0.  5 

1.0 

6.0 

WT-VI 

Arapaho  C  with 

Small 

3.  425 

2.  3,  4,  5 

1.0 

1.0 

5.  0 

nose  cone 

Large 

3  425 

2,  3,  4,  5 

1 . 0 

1  0 

5.  0 

Blunted  elliptical 

Small 

3.  425 

2,  3,  4,  5 

0.  5 

0.  5 

5.  0 

cone 

Large 

3.  425 

2.  3,  4.  5 

5 

0  5 

5  0 

Max 


8 

8 

9 

8 

8 

8 

8 

8 

8 

8 

10 

10 

7 

7 

7 

7 


TABLE  V  -  FOREBODY  P./P  RATIOS 

b  oo 


P  / 
b' 

P 

oo 

Forebody 

M  =2 
oo 

M  =3 
oo 

M  =4 

oo 

M  ^  5 

CD 

Baste  Arapaho  C 

0  56 

1 

0.  25 

0.23 

0.  24 

Arapaho  C  with  nose  cone 

o 

o 

0.24 

0.  19 

0.  17 

Blunted  elliptical  cone 

0.  49  I 

i _ 1 

0.  29 

i _ 

0.  17 

0.  15 
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Figure  JO  -  Drag  Coefficient  vertus  Viach  Number,  Ara.pa.ho  C  with 
Noee  Cone  Forebody,  Ac/A^  =  3.425 
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4.  REYNOLDS  NUMBER  EFFECTS 

In  all  cases  during  the  parachute  performance  tests,  the  Reynolds  num¬ 
bers  were  such  that  the  boundary  layer  at  the  base  of  the  forebody  was 
turbulent;  the  viscous  portion  of  the  wake  was,  therefore,  also  turbulent. 
Any  variations  in  parachute  drag  due  to  operation  in  a  laminar  viscous 
wake  could  not,  therefore,  be  determined. 

The  variation  of  parachute  drag  behind  a  given  forebody  depends  on  the 
profile  of  the  flow  field  entering  the  parachute  A  slight  variation  of 
velocity  profiles  with  a  variation  in  total  pressure  can  be  seen  in  Fig¬ 
ures  8  through  15.  For  the  wake  survey  tests,  the  temperatures  were 
held  constant  at  each  Mach  number;  the  variation  of  Reynolds  number 
can  be  directly  related  with  the  variation  of  total  pressuie.  Considering 
that  the  Reynolds  number  variation  for  the  parachute  tests  was  approxi¬ 
mately  half  the  variation  for  the  wake  survey  tests,  and  that  the  wake 
profiles  shown  in  Figures  8  through  15  do  not  indicate  significant  varia¬ 
tions  with  Reynolds  number,  it  can  be  safely  assumed  that  the  changes  in 
parachute  drag  over  the  range  of  Reynolds  numbers  tested  will  Le  very 
slight. 

From  the  w 4  *»d -tunnel  test  results  of  WT-HB,a  slight  correlation  between 
variations  of  Reynolds  number  and  variations  of  drag  coefficient  magni¬ 
tude  was  observed.  However,  inaccuracies  of  the  drag  measurement 
system  may  have  affected  the  repeatability  of  these  slight  'rends.  Fig¬ 
ures  32  through  35  illustrate  the  effect  of  the  free  stream  Reynolds 


Figure  32  -  Drag  Coefficient  versus  Unit  Reynolds  Number,  Basic 
Arapaho  C  Forebody  (M  =  2  and  3) 
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Figure  33  ~  Drag  Coefficient  versus  Unit  Reynolds  Number,  Basic 
Arapaho  C  Forebody  (M^  =  4  and  5) 
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Figure  35  -  Drag  Coefficient  versus  Unit  Reynolds  Number,  Blunted 
Elliptical  Cone  Forebody  (M  =  4  and  5) 

number  on  various  fovebody -parachute  combinations.  The  canopy  loca~ 
tion  of  X/D  a  6  was  chosen  because  of  the  near  optimum  performance 
exhibited  by  the  parachutes  at  this  location. 

The  figures  indicate  that,  behind  the  blunted  elliptical  cone  model,  there 
is  a  slight  increase  in  drag  coefficient  with  an  increase  of  Reynolds  num¬ 
ber.  However  the  Reynolds  number  trends  appear  to  be  reversed  for 
the  axisymmetric  forebody.  An  explanation  of  this  trend  reversal  has 
not  been  determined. 

5.  FOREBODY  EFFECTS 


Of  the  three  forebody  models  used  in  the  parachute  performance  testing, 
two  (the  b.  f  ic  Arapaho  O  and  the  Arapaho  C  with  nose  cone  forebodiesl 
were  axisymmetric;  the  third  (blunted  elliptical  cone)  had  an  elliptical 
cross  section  at  the  base.  The  third  model  was  tested  to  simulate  the 
wake  of  a  lifting  re-entry  body. 


The  forebody  effects  on  drag  coefficient  magnitudes  are  shown  in  Fig¬ 
ure  36.  The  data  shown  on  this  figure  were  obtained  from  the  second 
series  of  tests  (WT-YI)  and  reflect  a  high  degree  of  accuracy.  The 
data  indicate  the  effects  of  f'rebody  asymmetry  on  parachute  perfor¬ 
mance.  The  PARA5QNIC  performance  appears  to  have  a  higher  degree 
of  Mach  number  sensitivity  in  the  wake  of  the  asymmetric  forebody  as 
opposed  to  the  axisymmetr’c  forebody.  The  forebody  effect  is  most 
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apparent  at  Mq,  =  3.  At  Mqq  =  2,  the  effect  is  not  quite  as  apparent; 
however,  the  Mco  =  2  data  for  the  axisymmetric  forebody  may  be  in 
error  due  to  slight  damage  of  the  canopy  suspension  line  >.  This  damage 
is  discussed  further  in  Item  6.  below.  At  M^,  “  5,  the  low  drag  coeffi- 
cient  values  for  the  blunted  elliptical  cone  data  are  attributed  to  severe 
wake  modification.  This  also  is  discussed  further  in  Item  6. 

6.  CANOPY  AXIAL  LOCATION  EFFECTS 

The  axial  location  of  the  canopy  aft  of  the  forebody  was  found  to  become 
critical  at  the  extreme  ends  of  the  X/D  range  tested.  At  the  near  X/D 
locations  (X/D  =  4  to  4.  5),  the  proximity  of  the  canopy  to  the  forebody 
base  resulted  in  severe  wake  deformations.  In  some  cases,  the  modi¬ 
fication  was  of  such  severity  that  the  wnke  was  considered  to  be  "blown;" 
i.  e. ,  flow  leaving  the  forebody  bas''  separates  and  passes  around,  rather  than 
through,  the  canopv.  In  the  blown  wake  conditions  the  flow  field  entering 
the  canopy  is  subsonic  and  an  extremely  low  drag  coefficient  results. 

At  canopy  inlet  locations  greater  than  X/D  =  7,  the  oscillation  stability 
was  found  to  be  very  time  dependent  and,  at  times,  the  oscillations  be¬ 
came  divergent  in  nature.  The  breathing  action  of  the  parachute  also 
appeared  to  be  X/D  dependent,  with  the  magnitude  of  the  breathing  in¬ 
creasing  with  decreasing  distance  between  canopy  and  forebody  base. 

From  a  stability  standpoint,  therefore,  the  optimum  location  of  the 
canopy  inlet  was  found  to  be  at  approximately  X/D  =  6. 

Figure  30  indicates  that  at  the  X/D  -  6  canopy  inlet  position,  the  drag 
coefficient  is  also  optimized.  When  the  Mach  number  2  data  (which  do 
not  follow  this  trend)  are  examined  in  more  detail,  and  specifically  the 
oscillation  angles  associated  with  these  data  points,  it  can  be  shown  that 
the  upper  two  data  points  exhibited  oscillation  angles  varying  from  10 deg 
to  13  deg;  the  middle  two  data  points,  angles  from  6  deg  to  8  deg;  and 
the  lower  two  data  points,  no  oscillations.  The  three  data  points  for  the 
small-mesh  parachutes  at  the  Mach  2  conditions  were  affected  by  a  near 
failure  of  the  parachute.  Immediately  after  completion  of  the  Mach  2 
points,  inspection  of  the  small -mesh  parachute  indicated  some  slippage 
of  the  suspension  lines  at  the  point  of  attachment  to  the  swivel.  The  slip¬ 
page  of  the  suspension  lines  was  such  that  the  parachute  inletwas  no  longer 
normal  to  the  flow. 

The  drag  coefficient  data  presented  in  Figure  35  for  the  blunted  elliptical 
cone  model  also  indicates  a  near  opti.  um  canopy  location  at  X/D  =  6 
except  for  the  Mach  number  5  conditions.  For  X/D's  =  5  and  6  at  = 

5,  and  X/D  =  5  at  Mqc  =  4,  the  modification  of  the  wake  was  such  that 
the  parachute  inflated  with  an  elliptical  shape  normal  to  the  flow  with  the 
major  axis  of  the  elliptical  canopy  in  the  same  plane  as  the  major  axis 
of  the  elliptical  <orebody  base.  The  degree  of  ellipticity  was  most  severe 
at  the  Mqj  =  5,  X/D  =  5  test  point.  At  these  test  points,  the  canopy  was 
virtually  motionless. 

7.  PARACHUTE  TO  FOREBODY  BASE  AREA  RATIO  EFFECTS 

The  parachute  to  forebody  base  area  ratios  (A  /Ajj)  varied  only  during 
the  first  series  of  parachute  performance  test§  (WT  -UB),  and  the  data 
accuracy  may  have  influenced  any  slight  trends.  The  A  /A^  effects  can 
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be  eeen  on  Figure*  32  to  35,  It  doe*  appear  that  the  smaller  A  /A^ 
ratio*  exhibit  slightly  smaller  drag  coefficient  values  for  the  bfhnied 
elliptical  cone  data.  For  the  smaller  area  ratio  (1.99),  the  viscous 
portion  of  the  wake  iff  considered  to  have  a  noticeable  effect  on  the  drag 
coefficient  value.  The  trend  of  smaller  drag  coefficient  values  with  the 
smaller  area  ratio  is  not  apparent  in  the  Basic  Arapaho  C  data.  In  this 
case,  the  smaller  area  ratio  (3.425)  is  great  enough  that  the  viscous 
portion  of  the  wake  has  little  effect  on  the  drag  coefficient. 

The  primary  effect  of  A  /A^,  as  discussed  previously,  is  associated 
with  the  canopy  locationcat  which  wake  modification  is  encountered. 
Reference  32  indicates  that  the  X/D  position  of  the  parachute  behind  the 
forebody  at  which  wake  modification  occurs  is  directly  proportional  to 
the  ratio  of  the  parachute  size  to  the  forebody  size.  The  results  of  the 
wind-tunnel  tests  conducted  under  this  program  tend  to  substantiate  this 
statement-  At  the  Mqq  =  5,  X/D  =  6  conditions  (WT-UB)  for  the  blunted 
elliptical  cone  tests,  severe  canopy  breathing  and  wake  modification 
were  encountered  for  the  Ac/Ab  =  3.425  test  points.  However,  for 
A^/A^  =  1.99  at  the  same  Mach  number,  no  adverse  effects  were  en¬ 
countered  at  X/D  =  6,  and  only  slight  wake  modification  obser.ed  at 
X/D  =  5.  Similar  effects  were  observed  with  the  axisymmetric  fore- 
bodies;  however,  the  effects  were  not  nearly  so  severe. 

8.  GENERAL  STABILITY 

It  was  observed  that  with  the  decelerator  positioned  at  a  large  X/D  loca¬ 
tion,  the  oscillations  induced  in  the  parachute  became  divergent  in  nature, 
and  axial  repositioning  was  required. 

The  observations  led  to  the  conclusion  that  a  near  optimum  X/D  position 
should  exist  where  the  minimum  inflation -oscillation  instability  is  prr  - 
sent.  This  near  optimum  X/D  position  was  observed  to  be  at  X/D  =  6. 

The  drag  oscillation  frequency  (f^-)  does  not  appear  to  be  Mach -number 
dependent  within  the  range  tested.  The  canopy  oscillation  angles  (a)  and 
oscillation  frequencies  (f^)  show  a  decrease  in  magnitude  with  increasing 
Mach  number.  The  breathing  frequencies  (f^)  did  not  appear  to  be  Mach- 
number  dependent.  However,  behind  the  blunte  d  elliptical  cone,  the 
breathing  pulsations  were  time  dependent  and  usually  of  small  magnitude 
so  that  breathing  frequencies  could  not  be  determined.  This  indicates 
that  although  the  drag  coefficient  experiences  degradation  with  increas¬ 
ing  Mach  number,  the  stability  is  enhanced. 

The  canopy  oscillations  were  of  dif’erent  form  behind  the  two  forebodie». 
Behind  the  basic  Arapaho  C,  the  oscillations  were  predominately  of  the 
coning  type  in  which  any  point  the  riser  fine  follows  a  circular  path 
in  a  plane  normal  to  the  direction  of  flow.  Bshind  the  blunted  elliptical 
cone,  few  coning  oscillations  were  observed.  Oscillations  behind  the 
blunted  elliptical  cone  were  predominantly  in  the  vertical  plane.  How¬ 
ever,  in  some  cases  at  far  X/D  positions,  the  oscillations  .in  the  verti¬ 
cal  plane  were  combined  with  higher  frequency  oscillations  in  the  hori¬ 
zontal  plane  A  point  on  the  riser  line,  then,  traces  a  zig-zag  path. 
Unfortunately,  frequencies  and  amplitudes  cots l{i  not  be  determined  for 
oscillations  in  the  horizontal  plane. 
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Inflation  stability  for  parachutes  ■’  a  rating  behind  the  blunted  elliptical 
cone  was  highly  X/D  sensitive.  In  ^  t  all  cases  w**rh  the  large- 
diameter  parachute  positioned  at  X/  4.  5,  severe  canopy  breathing 
pulsations  and  wake  deformations  wt  r  neon  iered.  Behind  the  basic 
Arapaho  C,  at  the  same  X/D,  Ac/A^  r.  <  .me  with  the  same  roof  mesh, 
canopy  breathing  pulsations  and  wake  de.  ormities  were  less  severe. 

Throughout  the  tunnel  tests,  on-line  observations  indicated  that  the 
canopy  was  assuming  an  elliptical  shape  normal  to  the  direction  of  flow 
while  operating  behind  tne  blunted  elliptical  cone.  The  *  sree  of  ellip- 
ticity  increased  with  increasing  Mach  number  and  decre_  »ng  canopy 
X/D  location.  This  ellipticity  could  not,  however,  be  determined  from 
the  photographic  coverage. 

The  type  of  roof  mesh  was  found  to  have  little  effect  on  oscillation  am¬ 
plitudes  or  frequency,  but  to  have  an  effect  on  inflation  stability  (breath¬ 
ing  pulsations).  The  large -mesh  canopies,  being  much  more  *  sgid  tha»» 
the  small -mesh  canopies,  did  not  exhibit  the  high-frequency,  low-ampli¬ 
tude  canopy  flutter  observed  with  the  small -mesh  canopies.  The  large  - 
magnitude  canopy  pulsations,  referred  to  as  breathing,  were  encountered 
with  both  canopies. 


(Reverse  is  blank) 
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SECTION  V 
DRAG  PREDICTION 


1 .  GENERAL 

This  section  presents  a  method  for  calculating  the  pressure  drag  of 
PARASONIC  parachutes  operating  in  supersonic  axisymmetric  wakes. 
The  method  uses  the  local  wake  profiles  as  the  flow  environment  and  the 
flow  through  a  curved  shock  wav..  *o  determine  pressure  distribution  on 
the  parachute  canopy.  This  is  a  gx  'at  improvement  over  previous  meth¬ 
ods  (Reference  33),  which  assume  flow  conditions  at  the  parachute  inlet 
equal  to  the  free  stream  conditions  ahead  of  the  forebody  and  a  normal 
shock  standing  at  the  parachute  inlet.  The  canopy  flow  field  in  this 
analysis  is  determined  using  the  following  basic  assumptions: 

1.  A  detached,  curved,  »hock  is  present  ahead  of  the 
canopy  inlet. 

2.  The  presence  of  the  canopy  in  the  wake  produces 
negligible  distortions  of  the  wake  flow  field. 

3.  The  flow  is  sonic  through  the  roof  grid. 

4.  The  flow  stagnates  at  solid  portions  of  the  canopy 
roof. 

5.  Secondary  effects  due  to  separated  flow  and  skin 
friction  can  be  largely  ignored. 

2.  FLOW  FIELD  ANALYSIS 

The  flow  analysis  begins  with  the  determination  of  the  local  liow  para¬ 
meter  profiles  in  the  wake.  The  most  important  parameters  required 
for  this  analysis  are  the  local  Mach  number  (MjJ,  static  pressure  (Pj^), 
pitot  pressure  (Pt,  ).  and  dynamic  pressure  (qj*).  These  profiles, 
obtained  either  frofn  experimental  tests  or  suitable  analytical  prediction 
as  described  in  Section  III  of  this  report,  represent  wake  profiles  that 
are  free  from  any  adverse  effects  caused  by  the  presence  of  a  trailing 
body  in  the  wake  flow  field.  The  "free"  wake  profiles,  therefore,  cannot 
be  considered  valid  for  the  cases  in  which  the  trailing  parachutes  have 
been  known  to  modify  the  wake  ct’-ucture  radically.  Figures  37  and  38, 
traced  directly  from  schlieren  photographs  obtained  during  wind-tunnel 
tests  (Section  IV),  illustrate  the  effect  of  the  parachute  axial  positioning 
on  the  wake  flow  structure.  The  flow -field  modification  is  quite  evident 
with  the  parachute  inlet  positioned  at  X/D  *  4.  5.  With  the  parachute  in¬ 
let  positioned  at  X/D  =  8,  however,  the  wake  trailing  shock  does  not 
appear  to  be  deformed.  The  position  of  the  trailing  shock  indicates, 
however,  that  the  presence  of  the  parachute  causes  a  slight  widening  of 
the  wake  neck.  This  slight  modification  of  the  wake  is  not  considered 
in  the  drag  prediction  analysis 
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On  the  assumption  that  the  wake  profiles  have  been  determined,  the 
next  step  in  the  flow  analysis  is  to  select  profiles  at  an  axial  station 
immediately  ahead  of  the  canopy  oblique  shock.  When  this  procedure 
is  used,  the  influence  of  the  parachute  confluence  point  and  suspension¬ 
lines  on  the  flow  field  is  neglected.  This  simplifies  the  analysis  but 
may  influence  the  accuracy  of  the  predictive  method. 

Analyses  of  schlieren  photos  show  that  two  general  types  of  shocks  ap¬ 
pear  around  the  supersonic  parachute:  detached  inline  shocks  and  at¬ 
tached  shocks  at  the  canopy  inlet.  The  shape  and  stand-off  distance  of 
the  inline  shock  is  determined  by  the  local  Mach  number  profile  and  the 
aerodynamic  porosity  of  the  canopy  (Reference  34).  The  shape  of  the 
attached  shock  is  determined  by  the  shape  of  the  canopy  and  by  the  Mach 
number  proiil:  behind  the  parachute  inline  shock  and  is  p-  esent  only 
when  the  flow  at  the  inlet  edge  is  supersonic.  The  shap  of  the  attached 
inlet  shock  need  not  be  known  for  the  drag  prediction  aria',  sis.  The 
shape  of  the  inline  shock  must,  however,  be  known  in  order  to  determine 
the  pressure  distribution  inside  the  canopy.  The  shape  of  the  inline 
shod:,  unfortunately,  is  not  easily  predicted.  The  shock  wave  oscilla¬ 
tions  encountered  with  relatively  low -porosity  canopies  (Reference  34) 
along  with  the  interaction  of  the  shock  wave  and  the  flexible  suspension 
lines  result  in  a  very  unstable  inline  shock  and,  therefore,  a  nonsteady 
flow  field  entering  the  canopy  inlet.  Inflation  and  oscillation  stability 
also  affect  the  stability  of  the  inline  shock.  The  flow  analyst*,  then, 
assumes  the  canopy  to  be  motionless  and  the  inline  shock  to  be  quasi - 
steady. 

Analysis  of  the  schlieren  photos  indicate  one  other  interesting  phenomenon 
associated  with  the  cai opy  inline  shock.  In  some  cases  a  seemingly  in¬ 
tense  and  steady  inline  shock  would  lose  its  identity  at  the  very  center. 
This  "hole"  in  the  center  of  the  shock  indicates  subsonic  flow  conditions 
at  the  wake  centerline.  The  subsonic  core  at  the  wake  centerline  is 
most  probably  caused  by  the  formation  of  a  wake  by  the  parachute  attach¬ 
ment  swivel  and  confluence  ring.  In  cases  where  the  forebody  wake 
centerline  velocity  is  at  low  supersonic  velocities,  further  deceleration 
of  the  flow  by  the  swivel  and  confluence  ring  then  results  in  a  subsonic 
core  of  air  entering  the  parachute  inlet. 

3.  CANOPY  PRESSURE  E-iSTRiBUTiON 

The  pressure  distribution  on  the  supersonic  parachute  canopy  is  deter¬ 
mined  by  two  flow  piths:  internal,  through  the  canopy,  and  external, 
around  the  canopy.  The  fiow  field  in  which  the  canopy  is  engulfed  is 
influenced  bv  the  detached  oblique  shock  ahead  of  the  canopy.  With  free 
stream  profiles  obtained  from  experiment  or  from  the  method  described 
in  Section  III  and  a  detached  shock  wave  shape  determined  from  either 
experiment  or  some  applicable  predictive  methc-d.  fiow  parameters  aft 
of  the  detached  shock  wave  may  be  determined  by  the  use  of  oblique 
shock  relations.  Both  the  total  and  static  pressures  arp  required  to 
determine  the  pressure  distribution  on  the  canopy.  In  the  absence  of 
fu.ther  perturbations  aft  of  the  shock  waves,  the  parameters  are  com¬ 
bined  within  stream  layers.  If  the  shock  is  discontinuous  at  the  wake 
centerhne  region,  the  fiow  is  still  considered  isenlropic  and  is  treated 
as  one -dimensional  and  compressible. 
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The  axial  velocities  in  the  stream  layers  are  assumed  to  stagnate  on 
the  solid  portions  of  the  canopy  mesh  and  on  the  canopy  roof  cap  (see 
Figure  39).  The  surfaces  of  the  canopy  skirt*  however,  are  assumed 
to  be  acted  upon  by  the  static  pressure  of  the  stream  layer  entering  the 
parachute  at  the  inlet  edge.  The  parachute  canopy  acts  as  a  supersonic 
diffuser  in  that  the  flow  inside  is  decelerated,  and  the  kinetic  energy  of 
the  flow  is  converted  to  a  pressure  rise.  At  each  mesh  opening  the  flow 
is  assumed  to  be  sonic;  the  static  pressure  in  the  opening  is  then 

P  =  0. 528  P  . 
e  le 

The  external  flow  influence  on  the  canopy  is  divided  into  two  regions. 

The  first  region  is  skirt  surface,  beginning  at  the  inlet  edge  and  ending 
at  the  juncture  of  the  skirt  material  and  the  roof  mesh.  From  the  lead- 
ing  edge  of  the  skirt  to  the  point  of  maximum  camber,  the  skirl  is  as¬ 
sumed  to  have  a  constant  slope,  and  the  pressures  on  this  surface  are 
assumed  to  follow  conical  flow  theory.  Along  the  skirt  surface  aft  of 
the  point  of  maximum  camber,  the  flow  expands  and  is  assumed  to  he 
analogous  to  the  flow  over  a  rirg-type  airfoil,  whereon  the  pressure 
coefficient  along  the  upper  surface  may  determined  from  Reference  35  ar 

S  -  ?  •  V  • 

4.  PRESSURE  DRAG  COEFFICIENT  EVALUATION 

The  pressure  drag  associated  with  a  cupencmc  parachute  is  the  major 
drag  component,  and  its  analytically  determined  magnitude  is  an  accu¬ 
rate  estimation  of  the  total  drag  of  the  parachute.  Determination  of 
secondary  drag  components  is  beyond  the  scope  of  this  analysis.  The 
pressure  drag  is  evaluated  using  the  method  described  below. 

The  pressure  drag  is  evaluated  by  a  summation  of  the  axial  components 
of  the  pressure  differential  acting  on  the  parachute  canopy  (see  Figure  40). 

Examining  a  ring  ot  radius  r.  and  width  ds,  the  tangent  at  point  A  inter¬ 
sects  the  parachute  centerline  at  an  angle  d-  The  surface  area  of  the 
ring  is  then 


dA  ,  =  2»rds  , 
surf 


but 


therefore. 


d#  =  dr  esc  c  , 


dA  -  s  2f  esc  d  rdr 
auri 


The  pressure  d’ffc rertiax  acting  cj>n  the  ring  is 

2  *°int  '  • 
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For  regions  of  constant  pressure  differential,  the  incremental  pressure 
drag  coefficient  can  be  expressed  as 


AC 


DP 


AP 

^oo 


mm 


To  simplify  the  calculation  procedure,  this  equation  may  be  rewritten  as 


^'“'DP  internal  ~P  external' 


or 


iCDI'  =  ACP 


ft)  -ft) 

The  total  pressure -drag  coefficient  is  then 

n 

V  AC. 


'DP 


.mmmd. 

1=1 


RESULTS 


Using  the  methods  described  in  Item  4  above,  several  sample  calcula¬ 
tions  were  made  and  compared  to  drag  coefficient  values  obtained  in  ex- 
perimenta’  tests  (Appendix  V).  For  the  sample  calculations,  the  para  • 
chute  canopy  was  broken  up  into  11  annular  rings,  each  assumed  to  have 
a  uniform  pressure  differential  across  its  surface.  The  location  and 
size  of  these  integration  strips  was  based  on  the  construction  of  the 
parachute  models  tested  and  on  the  pressure  distribution  on  the  canopy. 
Figure  39  shows  the  canopy  construction  and  the  coordinates  of  the 
integration  strips. 

Strips  1  through  3  make  up  the  total  skirt  surface.  Strips  4  through  9 
comprise  the  mesh  portion  of  the  canopy  roof.  These  strips  are  divided 
into  open  and  closed  areas.  Strips  4  and  6  are  the  open  strips  and  equal 
the  total  open  area  of  the  canopy.  Strips  10  and  11  correspond  to  the  non- 
porous  cap  on  the  canopy  roof.  The  choice  of  1  1  integration  strips  was 
deemed  sufficient  for  obtaining  reasonable  results  from  the  initial  calcu¬ 
lations  . 

The  waj  ;  profiles  used  for  the  initial  calculations  were  obtained  from 
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wind-tunnel  tests  (Appendix  V).  The  wake  profiles  were  obtained  in  both 
horizontal  and  vertical  planes  during  the  wind-tunnel  tests.  However 
the  forebody  mounting  strut  was  in  the  horizontal  plane,  and  the  wake 
n  ..aimrements  in  that  plane  were  considered  to  be  influenced  by  the  wake 
of  the  strut.  The  wake  profiles  in  the  vertical  plane  were  then  used  and 
were  assumed  to  be  axisymmetric.  The  shock  shapes  were  obtained 
directly  from  schlicren  photographs  of  the  test  items.  From  the  shock 
shape  and  wake  profiles,  the  conditions  aft  of  the  shock  were  determined. 

A  first  assumption  was  that  the  streamlines  aft  of  the  shock  remain 
parallel  to  the  wake  axis.  The  results  of  computations  made  using  this 
assumption  indicated  that  accounting  for  any  radial  transposition  of  the 
streamlines  would  not  substantially  increase  the  accuracy  of  the  calcu¬ 
lation  method. 

Table  VI  lists  the  conditions  and  calculated  values  for  one  of  the  test 
cases.  In  Figure  41,  results  obtained  using  the  drag  predictive  method 
are  compared  with  experimental  wind-tunnel  data  (Section  IV).  Table  VII 
presents  the  same  data  in  tabular  form.  The  selection  of  canopy  axial 
positioning  at  each  Mach  number  was  based  on  parachute  stability  and 
the  availability  of  a  well  defined  shock  shape  from  schlieren  photographs. 
As  can  be  seen,  there  is  good  correlation  between  the  calculated  and 
experimental  values.  In  most  of  the  test  cases  the  calculated  values 
have  r, mailer  magnitudes  than  the  corresponding  experimental  data.  The 
accuracy  of  the  method  does  not  appear  to  be  highly  Mach -number  de¬ 
pendent;  however,  as  the  Mach  number  increases,  the  possibility  of  wake 
modification  due  to  proximity  of  decelerator  and  forebody  becomes  greater. 


TABLE  VI  -  SAMPLE  CALCULATION  DATA 


Wake  profile  parameters 


Y/D 

pi 

(psf) 

M1 

e 

Mi 

sin  8 

P2 

(psf) 

Pt2 

(psf) 

Cp2 

Cpt2 

C 

pe 

0 

28.40 

0.  70 

*  »  • 

e  •  » 

•  •  • 

30.4 

t  •  « 

.  0756 

»  •  • 

0-2 

28.  40 

1. 70 

44° 

1.  18 

41.4 

66.9 

.  0904 

.267 

.  048 

0.  3 

27.60 

1.92 

40° 

1.23 

44.  U 

68.8 

.  114 

.286 

.  055 

0.5 

28.80 

2.  70 

38° 

1.66 

87.8 

00 

.410 

.610 

.  228 

0.6 

28.  20 

2.76 

38° 

1.  70 

90.4 

119.01 

.431 

.630 

.  241 

_ I  111 _ J _ _ L _ I _ _ L__ _ _ 

♦These  data  were  calculated  at  the  following  free  stream  conditions:  (1)  para¬ 
chute  -  Mqq  =  3.  00,  qQQ  =  0.  98  psi,  Re/in  =  .  069  X  10  ,  X/D  =  6.  0,  = 

5.7  nai:  and  (2)  wake  measurements  -  =  3.00,  =  1.01  psi,  Re/in.  = 


W  • 

5.  7  psi;  and  (21  wake  measurements 
.069  X  106,  X/D  =  3.  Pt00  =  5.9  psi. 

+  Conical  flow  theory 

Ring  airfoil  theory 
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TABLE  Vn  -  EXPERIMENTAL  AND  CALCULATED  DRAG 
COEFFICIENT  VALUES 


M 

oo 

X/D 

CDP 

<cD) 

exp 

S/V  *  100 

exo 

2 

5 

0.  734 

mm 

86.4 

7 

0.873 

BBS 

92.  0 

3 

5 

0.  585 

0.494 

118.  3 

6 

0.  482 

0.  563 

85.  7 

4 

5 

0.  369 

0.423 

87.3 

7 

0.  361 

0.  340 

106.  0 

5 

6 

0.260 

0.  376 

69.2 

7 

0.  320 

0.  342 

93.  5 
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The  pressure  drag  is  the  major  drag  component, and  the  theoretically  de¬ 
termined  drag  magnitude  is  an  excellent  indication  of  the  supersonic 
parachute  performance.  Evaluation  of  other  drag -influencing  parameters 
listed  below  could  possibly  enhance  the  accuracy  of  the  predictive  method. 

1.  Skin  friction 

2.  Separated  flow  effects 

3.  Confluence  point  and  suspension  line  effects 

4.  Effective  porosity  variation  due  to  flow  compressibility 

Calculation  of  drag  components  attributed  to  these  phenomena  does  not 
appear  practical  from  an  engineering  standpoint;  therefore,  only  trends 
are  discussed. 

The  skin -friction  drag  can  be  neglected  without  sacrificing  appreciable 
accuracy.  The  drag  component  attributed  to  separated  flow  effects  may 
have  a  slight  effect  on  the  accuracy;  however,  evaluation  of  this  com¬ 
ponent  does  not  appear  practical.  The  influence  of  the  confluence  point 
and  suspension  lines  has  been  studied  (Reference  34),  and  results  indi¬ 
cate  that  the  number  of  suspension  lines  does  affect  the  internal  pressure 
distribution  on  the  canopy-  The  effect  of  the  confluence  point  and  sus  - 
pension  lines  has  been  taken  into  account  somewhat  in  that  the  shock 
shapes  used  for  the  initial  calculations  were  obtained  in  the  presence  of 
the  confluence  point  and  suspension  lines.  The  variation  of  effective 
porosity  with  Mach  number  does  not  appear  to  be  an  appreciable  lactor 
when  dealing  with  the  low  porosities  associated  with  supersonic  para¬ 
chutes. 

Recommendations  for  future  work  are  included  in  Section  VII. 


7* 


AFFDL-TR  -67-192 
Volume  II 


SECTION  VI 
CONCLUSIONS 


1.  A  method  of  solution  was  developed  to  predict  the  flow  field  properties 
as  well  as  the  boundary  layer  momentum  thickness  and  momentum  de- 
feet  over  axisymmetric,  asymmetric,  and  two  “dimensional  forebodies. 

2.  A  method  of  solution  was  developed  to  predict  the  properties  of  the  wake 
flow  field  behind  an  axisymmetric  forebody  in  a  supersonic  free  stream. 

3.  These  methods  have  been  programmed  in  FORTRAN  IV.  Results  using 
this  program  have  shown  good  agreement  with  experimental  wake  re¬ 
sults  over  a  large  range  of  free  -stream  Mach  numbers  and  axial  loca¬ 
tions  aft  of  the  forebody  base. 

4.  A  method  for  calculating  supersonic  par.  chute  drag  using  wake  profiles 
as  free -stream  conditions  was  developed. 

5.  Results  using  the  parachute  drag  predictive  method  have  been  compared 
with  experimental  data  and  good  agreement  has  been  shown. 

6.  The  results  of  the  PARASONIC  parachute  performance  tests  indicate 
that,  although  the  drag  coefficient  decreases  with  increasing  Mach  num¬ 
ber,  stability  generally  is  enhanced.  The  effects  of  wake  asymmetry  on 
parachute  performance  became  more  pronounced  as  the  Mach  number 
increased.  The  d.ag  coefficient  values  did  not  show  any  repeatable 
trends  associated  with  Reynolds  number  or  type  of  roof  mesh. 


(Reverse  is  blank) 
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SECTION  VII 
RECOMMENDATIONS 


The  following  task*  are  recommended  as  desirable  improvements  to  the  fore¬ 
body,  wake,  and  deceleraior  analyses: 

1.  Forebody  anaiyci* Improve  the  turbulent  boundary  layer 

momentum  defect  calculation  by  using  advanced  methods 

for  determining  body  friction  coefficients 

2.  Viscous  wake  analysis 

a.  Better  determine  the  starting  point  for  the  wake  calcu¬ 
lation  and  consider  modification  of  the  starting  point 
from  the  rear  stagnation  point  to  the  viscous  wake 
neck 

b.  Extend  the  analysis  to  flows  with  Prandti  number  not 
equal  to  1 

c.  Include  the  effect  of  axial  pressure  gradient  (this 
would  increase  the  complexity  of  the  analysis) 

d.  Develop  a  solution  for  the  near  wake  region  (this  would 
require  a  major  effort). 

3.  Inviscid  wake  analysis 

a.  Include  the  effect  of  radial  pressure  gradients 

b-  Account  for  the  effect  of  embedded  oblique  shock  wakes 
(wake  recompression  shock  and  flare  shock)  on  the 
inviscid  wake  profiles 

4.  Wake  solution 

a.  Investigate  methods  for  matching  the  viscous  wake 
and  mviscid  wake  edge  conditions 

b.  Apply  the  wake  solution  to  additional  experimental 
data  to  determine  the  choice  of  viscosity  model  and 
viscosity  coefficient  that  best  correlates  theory  and 
experiment 

c.  Perform  a  sensitivity  analysis  to  determine  which 
modifications  to  the  analysis  method  will  improve 
the  solution  the  most 

5.  Drag  prediction 

a.  Extend  the  drag  predictive  analyses  to  apply  to  non- 
porous  declarators 

b.  Conduct  experimental  tests  on  rigid  decelerator 
models  of  simple  geometry  (for  example,  cores  and 
spheres)  and  correlate  with  analytical  results 

(Reverse  is  blank) 
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APPENDIX  I 


BOUNDARY  LAYER  THEORY 


LAMINAR  BOUNDARY  LAYER  THEORY 

For  the  methods  developed  by  Kemp,  Rose,  and  Detra  (Reference  1)  for 
calculating  boundary  layer  properties  and,  in  particular,  calculating 
momentum  thickness  (0),  Stetson  (Reference  6)  has  presented  an  em¬ 
pirical  correlation  of  the  results  of  the  numerical  solutions  for  9,  which 
is  given  by 


-^4(0.491)0  -  0.090/3°-4)(peWe/Pw^w)0-386 


o  u  r 
e  e 


Here, 


-r 


o  u  u  r  ds 


and  s  is  the  distance  from  the  stagnation  point  as  measured  along  the 
surface  of  the  body;  k  =  0  for  two-dimensional  flow,  and  k  =  1  for 
axisymmetric  flow.  For  nonaxisymmetric  flow,  one  of  these  two  cases 
must  be  used  as  a  suitable  approximation.  For  example,  a  flat  delta  of 
small  thickness  ratio  is  best  approximated  as  a  two-dimensional  body, 
while  an  elliptical  body  with  its  major  and  minor  axes  on  the  same  order 
of  magnitude  is  more  suitably  approximated  as  a  quasiaxisymmetric 
body  (k  =  1). 

In  the  cbove  equations,  the  subscript  e  denotes  local  conditions  at  the 
edge  of  the  boundary  layer,  while  w  denotes  local  gas  conditions  at  the 
wall.  The  effect  of  the  pressure  gradient  parameter,  is  small  in 
Equation  24.  This  can  be  seen  since  0  ranges  from  one  at  a  two-dimen¬ 
sional  stagnation  point  to  Zero  in  a  xero  pressure  gradient  region;  thus, 
a  typical  value  ior  represents  no  more  than  a  10  percent  correction. 

In  Equation  24,  represents  an  effective  surface  distance  that  considers 
the  prior  boundary  layer  history.  For  purposes  of  calculation,  can 
be  rewritten  as 


(  <•>  =  P, 


ooR.T  R  <  ’ 

1  w 


where 


/s/R 

5*  7*  (if  4) 

S  CD 
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Pa  is  the  forebody  nose  stagnation  pressure,  Vqq  the free-stream  velocity, 
and  R  the  base  radius  for  an  axisymmetric  body  or  an  effective  base 
radius  for  a  no naxi symmetric  body)  assuming  y/RT  ■  constant. 


'he  boundary  layer  momentum  thickness,  9,  can  be  rewritten  in  terms  of 
The  result  is 


Here,  0  can  be  evaluated  from 


0  s  2  ^  *  con,t  •  (29) 

Then  the  local  boundary  layer  momentum  is  known,  M  can  be  evaluated 
from 


vX 


'circumference 


2a. 

P  u  eds 
e  e 


(30) 


whore  ds  is  a  circumferential  increment  of  length.  From  Equation  28, 

P  u26  can  be  written  as 
8  8 

”•“*  ’  />_V_r\1/V«_V/2/.>/#_«_\0.  386  MW  \V 


(^)  fe)  ttffeS 


an 


Thus,  the  integral  in  Equation  30  esn  be  evaluated  readily  ones  the  local 
flow  properties  have  been  used  to  evaluate  £  ,  8/13;.;  and  ^cu|9-  For  a 
circular  body  and  by  evaluating  the  integral  at  the  base,. 


M  » 

For  a  body  with  elliptical  cross  section*  M  can  be  *xpi<e**e<§  as 

/**/2  }  a  2 

« * 4  /  v:  *  G.R  rt*  • 
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where,  because  of  the  symmetry,  it  is  only  necessary  to  integrate  over 
one  quadrant«and  where 


Here,  b  and  a  are  the  semiaxes  of  the  ellipse.  is  unity  for  a  cylin¬ 
drical  body. 

Once  the  inviscid  flow  properties  along  the  edge  of  the  body  boundary 
layer  are  known,  and  assuming  the  wall  temperature  is  known  or  can  be 
estimated,  then  the  boundary  layer  momentum  thickness  and  integrated 
momentum  defect  can  be  evaluated  using  Equations  26  and  30. 

2.  TURBULENT  BOUNDARY  LAYER  THEORY 

Unfortunately,  the  state  of  knowledge  regarding  turbulent  boundary 
layers  is  not  so  well  developed  as  in  the  laminar  case.  Thus,  although 
there  have  been  many  approximate  treatments  of  turbulent  flow,  most 
of  which  can  be  supported  to  some  extent  by  using  existing  experimental 
data,  there  are  no  "exact"  theoretical  solutions  (such  as  those  of  Kemp, 
Rose,  and  Detra)  for  laminar  flow. 

One  approximate  solution,  the  results  of  which  lend  themselves  to  the 
present  problem,  is  that  of  Reshotko  and  Tucker  (Reference  71.  This 
solution  uses  the  momentum  integral  and  moment  -of -momentum  equa¬ 
tions  as  simplified  using  Stewartscn's  transformation.  To  solve  these 
»*o  equations,  a  skin  friction  relation  must  be  used;  Reshotko  and  Tucker 
cho-e  the  Ludwieg -Tillman  relation  in  a  form  suitable  tor  -romprea.8 >hie 
flow  with  heat  transfer  through  applying  the  reference  enthalpy  concept. 
The  above  equations  were  simpimed  further  by  using  an  approximate 
shear  stress  distribution  and  the  power  law  velocity  profile. 

The  moment -of -momentum  equation  ;s  needed  to  account  for  pressure 
gradient  effects  on  the  boundary  layer  velocity  profile.  The  method  of 
Reshotko  and  Tucker,  as  applied  to  insulated  surfaces,  is  quite  well 
founded.  However,  ion  noninsulated  or  nonariiabatic  wall  cases  the 
method,  though  qualitatively  correct,  is  based  on  some  speculative 
assumption*.  It  is  still  anticipated,  however,  that  for  such  <**«•*.  the 
method  will-yield  reasonable  quantitative  results..  The  .method  certainly 
»3.  re-fifes entalive  of  the  best  that  can  be  done  within  the  present  status 
c-1  turbulent  flow  theory, 

•ScShbtko  and  Tucker  m  Reference  7  present  their  result*  for  momentum 
thickness  in  i-i  a  .ransformed  momentum  thickness,  8j„.  The 

transforn -•■d'pi<M«=nturn  thickness,  9tr,  can  be  related  to  im»e  actual 
thickness,  bv 
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Here,  a0  and  uc  are  the  velocity  of  scurxi  and  the  kinematic  viscosity 
evaluated  at  the  stagnation  conditions  r(  the  local  external  stream,  re¬ 
spectively.  The  exponent  B  can  be  exp  ued  approximately  as 


1.2  +  1 .  *8  f  'Y“'\  -  1  • 

'  o  / 


where  T  is  the  mviscid  stagnation  r  nparature.  The  reference  tempera¬ 
ture,  r°  ,  is  given  ms 

iCi 

r  .  =  S.  50  (T  i  T  }  +  0.22  (T  -  T  }  •  {33) 

rei  v.  e  aw  :  e 

where  T  >*  the  adiabatic  wall  or  recovery  temperature 


/  1/3  V 

e  (i  +  f~  *%>  j  * 


T  a  T 
aw 


The  distance,  t  . .  in  Equation  36  is  the  starting  point  of  the  calculation. 
4.1  the  bo'tndxry  layer  is  turbulent  over  the  a&tire  body,  then  sj  -  0,  and 
cah  uiation  starts  at  the  forward  stagnation  point.  However,  it  the 
boundary  layer  is  in.  tally  laminar,  then  the  calculation  starts  at  s- s1# 
the  transition  point.  At  that  point, 

0,’urfcuient)  -  ©(laminar?  »  \4r>) 

The  previously  noted  equation  for  a  compressible  turbulent  boundary 
layer  (Equation  36).  cau  be  recast  in  terms  of  a  turbulent  distance  pa¬ 
rameter,  £r  ,  where 
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APPEND’S  II 

FOREBODY  FLOW  FIELD  PROGRAM 


The  forebody  geometry  is  stored  as  a  table  to  allow  all  necessary  geometric 
properties  to  be  determined  by  simple  calculation  from  the  geometric  coor¬ 
dinates.  The  storage  system  defines  the  body  profile  as  a  series  of  points 
that  are  the  values  of  the  distance  r  from  the  forebody  axis  to  its  surface  and  as 
a  function  of  the  axial  distance  x  from  the  nose  of  the  body.  Thirteen  pro¬ 
files  are  stored  on  13  planes  formed  by  rotating  the  vertical  plane  containing 
the  forebody  axis  about  this  axis  90  deg  in  7.5-deg  increments.  Mirror 
symmetry  is  assumed  about  the  ve-t;cal  and  horizontal  planes. 

The  one  geometric  input  required  by  the  boundary  layer  program  is  the  value 
of  r/R,  where  R  is  average  base  radius.  The  value,  r,  is  fr«und  as  function 
of  u  along  any  one  4  by  linear  interpcJationbetween  tM  points  stored  in  the  table. 

For  the  pressure  distribution  subroutine,  the  angle  of  the  forebody  surface  in 
e  cb  4  plane  with  the  body  axis  must  be  known.  To  calculate  this, 


r . 


V  4 


-  arc  tan 


L±_L_£ 


i  +  1.  * 


i,  4 


(43) 


The  application  of  the  Newtonian  or  tangent-cone  method  for  estimating  the 
static  pressr  e  distribution  is  governed  by  the  following  logic.  The  free- 
system  M*ch  number  and  the  static  pressure  (M^  and  p,^)  are  inputs  to  the 
pressure  distribution  subroutine.  For  each  4  profile,  the  following  checks 
and  decisions  are  made. 


The  nose  angle  is  calculated  fro.n  the  stored  geometry, 


0n  ,  =  arc  tan 

0,  p  xlf  j  -  x. 


:o,  4 


(44) 


This  value  then  is  compared  with  a  stored  curve  of  9C  versus  M^,  where  0_ 
is  tl  •'  critical  cone  angle  between  conical  flow  with  an  attached  conical  shock 
and  flow  with  a  detached  shock. 


Case  A 


er  .  >  e 

0,  6  c 


If  00  ^  is  equal  to  or  larger  than  0^,  then  the  shock  wave  is  detached  from 

the  nose,  and  static  pressure  is  estimated  by  the  Newtonian  flow 

pe  =  ps  sin2  eQ  ^  . 

where  ps  is  the  nose  total  pressure.  Since  there  is  a  detaentd  shock,  ps  is 
the  free -stream  impact  pressure. 
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Case  B 


,  <  9 
0,  p  '  c 


If  8c,  is  less  than  6C,  then  an  attached  conical  shock  exists,  and  both  the 
static  and  nose  total  pressure  are  evaluated  from  the  conical  flow  subroutine. 
For  both  cases,  the  local  Mach  number  is  given  by 


M 


V 


I  v  ~  i 


v  -  i 


where  p-Q  is  the  local  total  pressure.  In  this  region,  pQ  =  ps. 


This  solution  is  carried  out  at  each  point  on  the  body  profile.  As  long  as 
8is  ^  +  the  local  total  pressure  is  equal  to  the  nose  total  pressure 

for  both  cases.  If  the  shock  wave  is  detached,  the  static  pressure  is  obtained 
by  the  Newtonian  method  until  0^  ^  <  8C.  after  which  the  conical  flow  eolation 
is  used. 


If,  however,  9^  ^  <  +  l  £  a  corner  shock  forms,  and  the  local  total 

pressure  changes.  To  evaluate  this  condition,  the  flow  properties  in  the 
regionxj,  „  i  to  are  uac-d,  and  the  angle  i)  =  8^  +  i  ^  ^  -  8^  ^  is  the 
conical  serniangle  of  interest.  As  before,  there  arc  two  possible  cases. 


For  Case  A  where  0  £  8_  for  Mo 

C  Cl 


ei  I 


1  r 


where  pQ.  ^  *  local  total  pressure  behind  a  detached  (normal)  shock.  This 
pressure  is  given  by 
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For  Case  B  where  8  is  less  than  0c  for  Me.  ,,  the  static  and  total  pressures 


—  '•  -  — . .  '  — c — _  i  ’  -  -  - - -  r*  — 

are  taken  from  the  cc  leal  "'low  routine  using  the  local  properties  at  point  i 

as  inputs.  Again,  the  Mach  number  is  given  in  either  case  by 
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APPENDIX  ID 

BOUNDARY  BAYER  TRANSITION  PREDICTION 


1.  GENERAL 

A  criteria  for  boundary  layer  transition  was  established  that  could  be 
incorporated  into  a  general  laminar -turbulent  boundary  layer  calcula¬ 
tion  program.  Any  such  criteria  must,  because  of  what  is  known  regard¬ 
ing  transition  and  turbulent  flow,  be  considered  as  leading  to  only  an 
approximate  estimate. 

;  i  developing  transition  prediction  criteria,  recognition  must  first  be 
made  of  the  factors  that  influence  transition,  the  parameters  in  terms 
of  which  transition  correlations  usually  are  made,  and  which  of  these 
factors  are  important  and  can  be  considered  realistically.  In  general, 
boundary  layer  transition  depends  on  the  Reynolds  number  of  the  flow, 
surface  roughness,  free -stream  turbulence,  noise,  compressibility  ef¬ 
fects  and/or  pressure  gradient  effects,  and  whether  the  boundary  layer 
is  being  heated  or  cooled.  Except  for  effects  due  to  surface  roughness 
and  free -stream  turbulence,  one  might  think  that  the  other  effects  are 
of  such  a  nature  that  methods  of  predicting  their  associated  influence  on 
transition  could  be  developed.  However,  the  effects  of  roughness  and 
turbulence  each  may  couple  with  the  other  factors,  thus  making  any 
attempt  to  separate  one  effect  from  another  extremely  difficult. 

As  an  example,  the  effect  of  surface  roughness  depends  on  the  thickness 
of  the  boundary  layer,  but  the  thickness  of  the  boundary  layer  depends 
on  the  compressibility  of  the  fluid  (as  manifested  through  the  local  Mach 
number)  and  whether  the  boundary  layer  is  being  heated  or  cooled,  as 
well  as  many  other  factors.  The  free -stream  turbulence  also  may  vary 
with  tunnel  conditions.  Thus,  the  individual  factors  do  couple  with  one 
another. 

In  evaluating  factors  affecting  transition  (see  Figure  42),  Low, in  Refer¬ 
ence  36,  presents  a  curve  showing  the  effect  of  free -stream  turbulence 
on  transition. 

Gazley  (Reference  37)  found  that  laminar  boundary  layer  oscillations  and 
their  damping  or  amplification  could  be  detected  only  when  the  free- 
stream  turbulence  was  less  than  0.001.  For  greater  turbulence  values, 
the  boundary  layer  oscillations  became  difficult  to  identify  because  of 
the  near  coincidence  of  theii  appearance  with  the  point  of  transition  to 
turbulence.  Gazley  also  reported  that  the  effects  of  free -stream  tur¬ 
bulence  tend  to  mask  out  the  effacts  of  other  variables .  This  may  ex¬ 
plain  the  disagreement  between  Stainback  (Reference  38)  and  Gazley. 
Staiaback's  results  tended  to  confirm  the  invariance  of  transition  Rey¬ 
nolds  number  with  Mach  number  up  to  Mach  8;  Gazley' s  results  indicated 
an  increase  of  transition  Reynolds  number  with  increasing  Mach  number. 

Both  Deem  (Reference  39)  and  Stewart  (Reference  40)  indicate  an  increase 
of  transition  Reynolds  number  with  an  increase  of  unit  Reynolds  number. 
Deem  also  observed  a  transition  reversal  phenomenon  due  to  changes  in 
wall  to  adiabatic  wall  temperature  ratios. 
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Figure  42  -  Free  Stream  Turbulence,  U'/U 

Correlations  cf  boundary  layer  transition  usually  are  presented  in  the 
form  of  a  transition  Reynolds  number  as  a  function  of  other  flow  simi¬ 
larity  parameters.  Although  experimental  data  frequently  are  presented 
using  the  free-stream  Reynolds  number  based  on  X  (Rex),  Reynolds 
numbers  based  on  momentum  thickness  and  displacement  thickness  also 
are  used.  If  the  transition  process  itself  is  viewed  as  originating 
through  local  boundary  layer  separation  occurring  because  of  a  local 
momentum  deficiency,  then  the  local  Reynolds  number  based  on  mo¬ 
mentum  thickness,  Reg,  would  seem  to  be  the  appropriate  dependent 
variable.  For  an  incompressible,  flat -plate  boundary  layer,  the  critical 
value  of  Reg  corresponding  to  the  boundary  layer  becoming  unstable  to 
small  disturbances  (e.  g.  ,  Tollmien-Schlicting  waves)  is  163  (Reference 
41).  On  the  other  hand,  the  transition  value  may  range  as  high  as  1300, 
depending  on  the  exact  surface  roughness  and  free-stream  turbulence, 
both  of  which  tend  to  decrease  the  value  of  the  transition  Reynolds  num¬ 
ber.  For  compressible  flows,  the  Mach  number  and  the  ratio  of  the  wall 
enthalpy  to  the  adiabatic  wall  enthalpy  also  become  important. 

Pressure  gradient  effects  (a  favorable  gradient  stabilizes  the  bound„ry 
layer  while  an  unfavorable  gradient  destabilizes  it),  although  important, 
are  much  more  difficult  to  consider.  The  unit  Reynolds  number  also  is 
included  in  many  correlations  of  transition  data.  Originally,  this  para¬ 
meter  mdy  ^ave  been  related  to  wind-tunnel  produced  disturbances. 
However,  there  is  now  some  evidence  that  it  also  is  an  important  effect 
in  correlating  ballistic  range  and  flight  data.  Its  role,  however,  is  not 
understood  from  a  physical  viewpoint. 
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2.  ANALYSIS 

In  developing  a  rather  simple  correlation  for  the  transition  Reynolds 
number  based  on  momentum  thickness,  only  the  effects  of  Mach  number 
and  hvv/haw  are  considered  directly.  Since  the  present  boundary  layer 
computer  program  calculated  Reg,  it  is  convenient  to  treat  the  transi¬ 
tion  problem  in  terms  of  this  parameter.  Furthermore,  the  free -stream 
Mach  number  range  of  interest  in  this  study  is  from  2  to  5.  For  a  blunt 
nose  body,  this  means  local  Mach  numbers  of  less  than  3.5. 

For  wind-tunnel  tests  on  Project  EUREKA,  the  wall  conditions  were 
such  that  hw  «  haw.  For  free -flight  calculations,  l^/h^  could  range  as 
low  as  0.  12,  based  on  an  upper  limit  of  6  for  Mqq.  Effects  due  to  surface 
roughness  and  free -stream  turbulence  will  not  be  included. 

In  general,  the  effect  of  the  local  Mach  number  on  the  transition  Rey¬ 
nolds  number  in  terms  of  X  can  be  represented  as  shown  in  the  following 
sketch  (Reference  42  and  43). 


A  curve  of  this  type  can  be  represented  by  a  power  series  of  the  form 

m 


Re 


■  r%L  5 


(47) 


If  only  the  first  four  terms  in  the  powe*-  series  are  used,  the  following 
boundary  conditions  can  be  applied; 


1.  M  =  0  Re,,  =  /ReY  \  (a  value  of  2.  5  X  10  was 
chosen)  *tr  \  tr/inc 

2.  M  =  0,  dReY  /dM  =  0 

Atr 

3.  M  =  M  , ,  dReY  /dM  =  0  (M  .  based  on  experimental 

data)  1  Xt-r  ‘ 


4.  M  = 

pe 


‘-Mi- -vTx) 

erimenta!  data)  ' 


(matched  with  ex- 
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If  boundary  conditions  3  and  4  are  selected  based  on  References  43  and 
44,  Equation  47  then  becomes: 


Re„  »  2.5  X  106(1  -  0.  12M2  ♦  0.  023  M3)  . 

*tr 

Since  Re  *  0.664  Re  (Reference  41),  then 

0  X 


(48) 


2  3  l/Z 

Re  w  1000  (1  -  0.  12  M  +  0.023  M  ) 
wtr 


(49) 


The  effect  of  h^Aaw  is  illustrated  in  the  following  sketch: 


This  curve  can  be  represented  by  a  power  series  of  the  form 

m 


Re 


tr 


(50) 


which  can  be  fitted  to  available  experimental  data.  When  the  range, 
0.2  .<  h  Amw  ^1.0,  i*  considered,  the  boundary  conditions  become 
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h  A 
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1.  Rev-  =  /Re Y  \  j  -  T 
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hj.  Re^y  (matched  with  experimental  data) 
h.,  dRev  /dh  »  0 
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0.  2,  ReY  (matched  with  experimental  data) 
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becomes 


and  where  Re^ 


“R'x 


the  correlation  of  Equation  49 
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Re.  *  1000  (1  -  0.  12  M  +  0.023  M  ) 

V 


1/2 


tr 


_  2  1/2 

(-2.29  +  17.  38  h  -  18  h  +  3.91  h  ) 


(51) 


where  h  =  Here,  the  effects  due  to  h  /h  are  based  on  the 

experimental  data  of  References  45  through  47. 

The  nature  of  this  correlation  is  that  transition  reversal  is  included. 
However,  the  second  transition  reversal  frequently  found  at  values  of 
1^/h  of  approximately  0.2  is  not  included.  The  correlation  (Equa¬ 
tion  51)  thus  is  limited  to  h  /h  >  0.2.  Reshotko  (Reference  48)  has 
qualitatively  explained  bothTransition  reversal  phenomena.  However, 
a  general  agreement  among  the  available  experimental  data  does  not 
exist;  the  correlation,  though  qualitatively  correct,  is  quantitatively 
only  an  approximation. 

In  representing  the  effects  of  Mach  number  and  h  /h  using  power 
series,  the  coefficients  in  th^se  series  should  be^lreaftd  as  functions 
of  surface  roughness  and  free -stream  turbulence. 

Pressure  gradient  efiects,  as  noted  before,  also  may  be  important.  On 
this  project,  the  primary  bodies  of  interest  are  such  that,  excluding  the 
nose  region,  the  body  is  divided  into  regions  of  zero  pressure  gradient. 
This  should  be  a  reasonable  approximation,  even  in  the  prediction  of 
transition,  with  the  one  exception  being  at  the  junction  of  the  cylindrical 
afterbody  and  the  flare  at  the  base  of  the  body.  The  adverse  pressure 
gradient  there  may  produce  transition  if  the  Reynolds  number  is  greater 
than  the  critical  value.  Thus,  it  was  suggested  that,  in  programming  the 
boundary  layer  calculations,  Reg  at  the  initiation  of  the  flare  be  checked 
against  <Re0)CRi'1'  a8  given  below: 


(He0) 


163  (1  -  0.  12  M  +  0. 023 


M3) 


1/2 


CRIT 


(-2.  29  +  17.  38  HT 


3.91  K3) 


1/2 


(52) 


This  equation  is  merely  Equation  51  with  the  incompressible  transition 
value  of  a  1000  replaced  by  critical  value  163).  If, at  the  flare,  Re0 
Re0  then  transition  is  assumed  to  If  Re8  <  ^e®CRlT*  t*'en 

the  nounoary  layer  remains  laminar.  Wit  i  the  inclusion  of  tm*  criterion 
in  the  program  boundary  layer  transition  should  at  least  be  accounted 
for  in  a  quantitatively  correct  manner. 

The  presebt  correlation.  Equation  52,  has  been  compared  with  experi¬ 
mental  results  from  Reference  40  on  the  effect  of  Tw/  Ta^,  obtained  at 
M(I)  -  3.  t.  Approximate  agreement  has  been  found  to  exist.  The  results 
of  Stetson  and  Rushton  (Reference  49)  also  have  been  compared;  these 
results  were  obtain*  1  blunt  cones  at  a  free -stream  Mach  number  of 
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Mqj  =  5.5  (see  Figure  43).  Here,  X-p  is  the  transition  distance,  X8W  is 
the  distance  for  the  high  entropy  .layer  to  be  swallowed,  ard  Reg-j.  is  the 
transition  Reynolds  number  based  on  momentum  thickness.  The1  pre¬ 
diction  based  on  the  present  method  is  shown  only  for  one  condition; 
however,  reasonable  agreement  exists. 

Equations  51  and  52  have  been  programmed  as  a  part  of  the  boundary 
layer  program  and  are  used  to  determine  the  approximate  location  of 
transition  if  the  exact  location  of  transition  is  not  specified  as  an  input. 
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APPENDIX  IV 
VISCOUS  WAKE  ANALYSIS 


1.  GENERAL 


An  integral  method  of  solution  was  applied  to  the  problem  of  determining 
the  properties  of  both  laminar  and  turbulent  viscous  wakes,  particularly 
to  the  problem  of  the  wake  flow  preceding  a  trailing  aerodynamic  de- 
celerator.  The  important  assumptions  contained  in  the  theory  are  as 
follows: 

1.  The  fluid  is  thermally  and  caloricaily  perfect, 
and  there  are  no  chemical  reactions  or  other  real 
gas  phenomena. 

2.  Prandtl  number  is  unity  (Pr  =  l). 

3.  Prandtl's  concept  of  viscous  flow  phenomena  is 
valid  for  the  high-speed  compressible  wake  in 
such  a  way  that  gradients  in  the  streamwise  or 
axial  direction  are  much  smaller  than  those  nor¬ 
mal  to  the  wake  axis,  and  boundary  layer  type 
equations  can  be  used. 

4.  The  details  of  the  base  flow'  and  free  shear  layer 
regions  can  be  largely  ignored.  Thus,  the  region 
of  ilidity  for  the  present  analysis  must  be  con¬ 
sidered  to  extend  from  somewhere  in  the  vicinity 
of  the  wake  neck  or  further  downstream  and  on 
down  into  the  far  wake. 

5.  Effects  due  to  the  cxistc-n  o  of  an  external  pressure 
gradient  are  negligible. 


By  using  these  five  assumptions,  the  momentum  integral  equation  can.  be 
written  as 


jB  j  U 


fh 

I  {U(u. 

Jq 


ulrdr  =  constant  with  respect  to  X 


(S3) 


Here,  0  is  the  momentum  thu  kness;  Equation  Si  simply  states  that  the 
wake  momentum  defect  is  constant  :r.  .‘he  axial  direction. 

By  ignoring  the  details  of  the  base  flow-  and  shear  layer  regions  and  as¬ 
suming  that  no  external  wor.».  ;s  done  or.  the  wake  fluid  m  these  regions, 
then  the  constant  n  Equation  S3  can  be  evaluated  by  equating  the  mo- 
rocaSUir.  defect  the  boundary  layer  at  the  base  of  tru  forebody  to  the 
momentum  defect  of  t nr  wa*r  can  be  stated  mathematically  as 

*‘bub^bd)  =  •  (S4) 
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where  the  subscript  1  refer*  to  properties  associated  with  the  outer  edge 
of  the  wake,  and  the  subscript  B  refers  to  properties  at  the  base  of  the 
vehicle,  where  the  boundary  layer  separates.  „  and  u  are  inviscid 
properties,  while  *>  is  the  forebody  base  diameter. 

The  continuity  and  momentum  equations  governing  this  Laminar  axisym- 
metric  high-speed  wake  problem  are 

’  0 

I# *  ,v  $7  *  7&!r"$7’ 

If'3-  l‘*< 

Here,  v  is  the  velocity  in  the  radial  direction,  r;  u  is  the  velocity  in  the 
axial  direction,  X.  For  laminar  flow,  the  properties  in  Equation  IS  are 
the  conventional  local  properties.  For  turbulent  flow  u,  v,  p,  and  p 
must  be  interpreted  as  time-averaged  quantities  (as  opposed  to  the  fluc¬ 
tuating  component  of  each  of  these  properties).  The  viscosity,  */, 
takes  a  different  form  for  laminar  flow  than  for  turbulent  flow.  The  axial 
direction  momentum  equation  can  be  expanded  and  expressed  as 

+  pvIt  =  +7  "-$7  •  <56) 

To  solve  this  equation,  the  Dorodnitzyn  transformation  is  introduced 

where 


ndn 


In  this  new  X,n  coordinatt  system,  the  transformed  wake  thickness  A 
is  found  through  integration  s~  that 


&  is  termed  the  wake  Thickness  bvt  is  really  the  wake  radius.  A  non- 
dimensional  variable,  1  -  n/A,  then  can  be  introduced  so  that  n  -  0  at 
the  center  of  the  wake  and  ?  =  l  at  the  outer  edge  where  Y  =  &  ,  Sr.  the 
X,  y  co-ordinate  system,  it  is  assumed  that  the  velocity  profile  car.  be  ex¬ 
pressed  as 
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'he  boundary  conditions  to  bo  used  then  are  a*  follows. 

i.  At  the  cento dine  of  the  wake,  Y  s  7  ~  0  and 


l.  At  the  centerline  of  the  wake,  Y  -  0  -  0  and 
Bn  _  jBu  __  - 

*37  "  "37  * 

j.  At  the  outei  edge  of  the  w;*ni,  Y  =  6  or  7  -  1, 

U  =  U  ,  . 

4,  At  the  outer  edge  of  the  wake,  Y  <5  or  ij  ~  5 

.  du  _  0 
~W~  Tv  '  0  • 

5,  At  the  centerline,  the  axial  direction  momentum 
equation  applies. 

6,  At  the  outer  edge  of  the  wake,  Y  -  6  or  n  - 
B^u  8*~u  n 

"  a7 

The  first  boundary  condition  yields 


Boundary  conditions  Z  through  4  yield 


a  j  -  0 


a  0  t  a  1  +  a  ?  +  a  x  ♦  a  d  -  l 


1  a ,  *  3  a  ,  *  4  \  -  0 


Boundary  condition  5  can  be  evaluate.!  cnl-  alter  the  a>.ia'  direct  101 
tnentufu  equation  -.s  t rausiorr.  e-i,  This  leads  t o 


„  T  *■  *•*  t  ■% 

C5  U  ^  ■  \  v  ^  y  : 

3x  15  *]£?*  7-  ~  o 


!.  t  n 


r><  cause  sv{  sr?  ayrunsetr  » 
recti  on  n_e  zero  On  '  1 :  ■ '  - 


m  tr.e  rauiai 


.  jis^ — 3Sfc  te*-*-*^ 
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leads  to 


2  pu  2 

“  ’  hr 


Evaluating  the  above  momentum  equation  at  the  wake  centerune  results  in 

du„  fi,  „  /  „  \  f>,  , 


„  0  _  u  1  p  8u\  ,  L  JlUm\  i0u 

0  dX  1  1  Vv  ;  r«  3n 


or 


du.  P,  ,  .  [A  2  r,  , 

„  a  — a  =  tr,  U_£  +  /I  Ou„ 

0  dX  1 


*  f  \  v  a  ,  <  *■  uu , 


Substituting  iro:n  Equation  59  results  in 


ruQ  _  Vl  . 

-Vo  dX  "  2  4a2  ' 

lA 


(61) 


Thus , 


pi\duo 


4a  =  l£/A2)  (ILO. 

2  vA)W,^ 


and  £' ^  oar,  he  expressed  as 


where 


(62) 


/  V.2.  /P1\H-U0 

>&)■<* 


(63 1 


The  first  fi .  e  ^boundary  conditions  thus  vie  id  tile  following  values  for  the 
coefficients  a^,  aj,  „  .  . 


ftn  ~ 


£ 

4 


!  00 
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The  velocity  profile  in  the  *ransformed  plane  then  becomes 
u  ”  S  P  2  '<  4  uo  a  4 

— r-1  =  -  2 r  +  O  +  (1  -  [in5  -  3 n  ) 

‘i  -  ui 

When  boundary  condition  6  is  used,  it  can  be  shown  that 


A  ■  (■  If) . 


The  velocity  profile  then  becomes 


u  ,/  2  Q  3  4  , .  £ 

=  (pr/  -8?  +  3t?  -  1}  -j- 


■  — +  i 
<  r  • 


For  £  =  24,  then  ~  0  and  X  =  0;  for  p  -  0,  then  =  1  and  X  =  oo. 

With  the  velocity  profile  now  known  as  a  function  of  P,  the  next  step  is  to 
relate  £  and  X.  This  can  be  done  by  using  Equation  63.  However,  to  do 
this  a  form  for  the  viscosity  coefficient,  //q,  must  be  specified.  Its  form 
will  depend  on  whether  the  flow  is  laminar  cr  turbulent;  the  evaluation  of 
the  relationship  between  p  and  X  will  be  discussed  in  the  next  section, 
where  the  various  special  cases  are  treated. 

ft  is  of  interest  though,  before  turning  to  the  solution  of  the  energy  equa¬ 
tion;  to  consider  the  relationship  between  G  and  A.  This  can  be  done  by 
transforming  Equation  53  so  that 

r  ,66’ 

Substituting  for  —  in  terms  of  rj  and  £,  there  results  upon  integration 
U1 


§i  j_  nr  _  _±_  pz 

a2  120  "  60, 480  "  120  '  5500  * 
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Where  there  is  no  pressure  gradient  in  the  wake,  0  is  constant  and  can 
he  evaluated  from  Equation  64;  then  Equation  6?  is  a  relationship  between 
A  and  £.  The  location  of  the  wake  neck  can  be  sought  by  differentiating 
Equation  67  with  respect  to  X.  This  results  in 


,  82  dA  /  1  2£\d£ 

■2TTd5T'  It TJS’imiw  * 


At  the  w&ke  neck. 


dA 

dX 


=  0. 


However,  d£/dX  ^  0;  thus,  the  location  of  wake  neck  is  determined  from 
the  condition 


1 


or 


£=  22.9. 


5500  120 

At  this  location,  the  centerline  velocity  then  will  be 


(68) 


u 


=  0.045. 


1 


(69) 


and  the  wake  thickness  or  radius  is 


(f)  =  3’24’ 

V0/neck 


(70) 


Although  the  above  results  are  of  interest,  they  are  only  approximate  since 
in  the  neck  egion  thp  effects  of  pressure  gradient,  in  both  the  normal 
a:-’'*  axial  directions,  may  be  important. 

In  ti  e  energy  equation,  the  high-speed  compressible  wake  can  be  written 
as 


9h  ,  ah 

/,u"3x  *  37 


r  -w— |  +  i 


c 

P 


(71) 


Multiplying  the  axial  direction  momentum  equation  by  u  and  adding  it  to 
the  above  results  in 


'$}  *  v£[(>  -  iM- $  •  <"» 


In  either  one  of  these  equations,  the  properties  must  be  considered  as 
local  fluid  properties  for  laminar  flow  and  as  time-averaged  properties 
tor  turbulent  flew. 
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For  Pr  =  1,  a  solution  to  the  above  energy  equation  in  terms  of  total 
enthalpy  is 


H  =  CjU  +  c^, 

where  the  constants  c.  and  c_  are  determined  from  the  boundary  or.  idi- 
tions  1  L 


u  =  u^ ,  H  =  H j 


u  =  u0.  H  =  H0 


The  total  enthalpy  profile  thus  can  be  expressed  as 


u.  -  u 

H,  -  H  =  (H.  -  HJ  -i— _ 

1  1  0  U;  -  U, 


or 


Hj  -  H 


i  - 


u 

U, 


1  - 


(73) 


(74) 


Hq  can  be  determined  from  the  condition  that  the  total  er  halpy  defect  is 
constant  in  the  wake  and  is  equal  to  that  of  the  vehicle  boundary  layer  at 
its  base.  For  the  wake,  the  total  enthalpy  defect  is 


r 

lUlHj  =  ZrT J0^u(Hl  *  H)rdr 


or  transforming 


(75) 


(76) 


Substituting  from  Equation  7  4  and  integrating  leads  to 
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I 


or 


1  - 


H, 


u. 


e 


2 

l 


(77) 


t 


i 


j 

i 

j 


f 

) 

s 


i 


f 

i- 

I 


For  the  vehicle  boundary  layer,  the  total  enthalpy  defect  is 


p  u  F 

S  3 


and  for  Pr  =  1,  this  can  be  reduced  to 


(78) 


(79) 


where  h  is  based  on  the  forebody  skin  temperature  at  the  base.  If  the 
total  enthalpy  defect  of  the  wake  is  equated  to  that  of  the  vehicle  boundary 
layer  at  the  base, 


PlUjKj  =  psUsHs(’7D0Tg)  • 

If  Equations  54,  77,  and  79  are  substituted,  then  it  follows  that 


(80) 


Thus,  knowing  the  velocity  profile  (Equation  65),  the  total  enthalpy  pro¬ 
file  can  be  determined  using  Equations  74,  77,  and  80.  The  static  en¬ 
thalpy  profile  readily  follows. 

To  transform  the  results  back  iniothe  physical  X,Y  coordinate  system,  use 
must  be  made  of  the  expression 


.104 
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4  =  2  -r/>dn  =  2  i 

&  J  n  *  n 


0  '0  1 

Substituting  for  the  enthalpy  profile  and  integrating  results  in 


=  2 


H1  '  H0 


■.&-?) 


/3  4  8  5  16 

l) 

t  2 

_£.  +  £_ 

[J'1  ‘S'7  +2n  - 

2/ 

24  2 

H, 


H1  'H0 


l1 


/  U_l 
2  '  2h, 


at  +  Alir,*  .*n5  +  £.£)]. 

2  12  U  5  c  2  /J 


(81) 


(82) 


2h, 


x 

24' 


(i- 


3,4+J|-^45,6.26,7tif/-^,9  + 

!>  toy 


10  / 


Thus,  except  for  the  relationship  between  £  and  X,  the  necessary  equa¬ 
tions  for  determining  the  properties  of  a  high-speed  compressible  wake 
have  been  derived  using  an  integral  technique  and  the  Dorodnitzyn  trans¬ 
formation.  For  a  given  X/D,  the  velocity  and  total  enthalpy  profiles  in 
the  transformed  coordinate  system  are  given  by  Equations  65  and  74  in 
terms  of  Ho,  the  centerline  total  enthalpy,  and  £,  a  form  parameter. 

Hq  can  be  determined  from  Equations  77  and  80  nnd  £  as  will  be  discussed 
in  the  next  section.  The  transformed  wake  thickness  or  radius  car.  be 
determined  from  Equation  67  and  the  physical  coordinate  r  from  Equation 
82.  To  obtain  numerical  results,  it  is  only  necessary  that  the  inviseki 
flow  field  be  specified,  and  that  an  estimate  be  made  of  the  vehicle  boun¬ 
dary  layer  momentum  thickness.  Although  some  cf  the  derived  equations 
may  seem  rather  cumbersome  in  terms  of  carrying  out  numerical  cal¬ 
culations,  the  results  are  of  an  algebraic  nature  and  in  principle  can  be 
readily  applied. 


It  should  be  re-emphasized,  however,  that  the  preceding  analysis  has  ig¬ 
nored  the  details  of  the  free  shear  layer  and  base  flow  regions.  The  re¬ 
sults  obtained  thus  must  be  considered  as  approximate;  their  range  of 
applicability  is  certainly  limited  to  the  region  downstream  of  the  wake 
neck. 


CONVENTIONAL  LAMINAR  WAKE 

To  evaluate  Equation  6?.  and  to-  obtain  £  as  a  function  of  X,  it  i?.  necessary 
to  specify  a  .relation  for  the  viscosity  coefficient.  For  laminar  flaw,  the 
viscosity  t  an  be  assumed  to  have  a  linear  dependence  on  the  fluid  tempera- 
tvire,  L  e. .  yu/RT  =  constant.  Since  in  this  anal vs is  the  radial  pressure 
gradient  is  aS9ctn..c  -to  be  zero,  -1  hen.  tfum  the  equation  of  state  it  e&n  be 
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shown  that  py  is  constant  with  respect  to  r  and  that/?^  -  Equa¬ 
tion  63  then  can  be  written  as  * 


"o  ,A2,  l‘o\  dU0 

'  “T*  \*i/"ar 


Since 


X  x,jL 

Uj  24 


%  !id i 

dX  “  24  dX  * 


Furthermore,  the  relationship  between  A^  and  £  is  known  from  Equa¬ 
tion  67.  Substituting  these  results  into  Equation  83  results  in 

tItT  =  -  (l  -  lr)^2T023x  •  (84) 

Since  P ^  has  not  been  specified  in  terms  of  £  or  X,  these  two  variables 
are  not  yet  uniquely  related  to  each  other.  To  do  this,  attention  must  be 
given  to  solving  the  energy  equation. 

Now  the  centerline  density,  can  be  expressed  as 

pl 

P°  ‘  RT0  "  Y-1  h0  ’ 

The  static  enthalpy  can  be  related  to  the  total  enthalpy  as 


h0  =  H0  -  2  ’ 


When  Equation  77  is  used,  pn  can  be  expressed  as 


#0  =  y  -  i 
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Since  the  relationship  between  u^/u.  and  £  is  known,  then  Equation  84 
takes  the  form 


120  5500  " 


(l  -  24)2* 


0lPl 


dX 


24/24  »l  v  -  1 


H1  - 


After  rearrangement,  this  takes  the  form 


Here,  X  is  measured  from  where  physically  the  wake  centerline  velocity 
is  zero,  i.  e.,  the  rear  stagnation  point.  For  known  flow  conditions, 
Equation  8f>  gives  the  variation  of  £  with  X.  When  the  results  of  Item  1 
are  used,  a  complete  solution  for  the  wake  flow  can  be  obtained. 

Equation  36  can  oe  recast  in  form  by  noting  that  the  bracketed  term  in 
the  denominator  can  be  written  as 


H. 


1  * 


and  that 


U1P1  1 _ 

A'i  Hj  y  -  1 


The  following  then  results; 
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where 


Re  -  V»° 
0  “  ^ 


and  is  based  on  wake  edge  conditions. 

Two  special  cases  are  of  interest.  One  is  that  of  an  adiabatic  forebody 
wall  such  that  8^  /of  -  0»  In  that  case,  Equation  87  becomes 


Re. 


Re, 


c 

hI 

(■4) 

(ii  £3 1 

\120  '  5500^ 

1  i 

(i  +  m2  —  X-i—L  m2  \ 

{  l  1  12  2  1  576/ 

d£-  (88) 


The  second  case  is  that  for  a  cold  A'all  where 


4 


l 


-  l. 


Here,  Equation  87  becomes 

rt 


Re, 


Re 


0 


J  2  4 


VI 20  5500/  \  2  1  24/ 


d£ 

24 


(89) 


Equations  88  and  89  indicate  that  for  £  small,  then 


X»£~  1 . 


Since  6  ~  A  for  £  small  and  similarly 


A2*£'1. 


Then, 


I 

6*X2  _  (90) 

£  small  corresponds  to  the  far  wake  region,  and  this  square  root  de¬ 
pendence  on  X  is  in  agreement  with  experimental  data  (Reference  13).  I* 
similarly  can  be  shown  that 


i  08 


1 


(91) 
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u 


U, 


and 

H0  -1 

1  -  rpocX  •  (92) 

1 

This  also  is  in  agreement  with  experimental  data. 

CYLINDRICAL  LAMINAR  WAKE 

The  nature  of  the  laminar  axisymmetric  wa  .e  in  the  presence  of  a  small 
diameter  wire  or  cable  at  its  centerline  is  such  that  the  wake  will  be  ap¬ 
proximately  cylindrical  (or  at  least  will  grow  very  slowly)  and  will  have  a 
zero  centerline  velocity  caused  by  the  no-slip  condition  at  the  cable. 

An  approximate  solution  corresponding  to  this  case  can  be  obtained  from 
the  results  of  Item  1  by  noting  that  Uq/Uj  equal  to  zero  is  equivalent  to 
£  equal  to  24.  Thus,  the  properties  of  a  cylindrical  wake  can  be  obtained 
by  setting  £  equal  to  24  in  the  previously  derived  equations.  £  WH1  not 
depend  on  X  for  this  simplified  case  since  there  will  be  complete  simi¬ 
larity  between  the  radial  profiles  at  all  X  direction  stations.  For  this 
case,  the  wake  velocity  profile  is  given  by 

~  =  *>nZ  -  8*3  +  3-7 4.  (93) 

U1 

and  the  transformed  wake  thickjiess.by 

£*  3.24.  (94) 

2 

where  0  is  determined  from  equating /qujGj  to  the  momentum  defect  of 
the  forebody  boundary  layer.  For  a  slender  body  with  a  zero  axial  pres¬ 
sure  gradient  in  the  wake  region,  9j  is  constant,  and  thus  A  also  is  a 
constant.  The  wake  thus  has  a  constant  diameter  and  is  cylindrically 
shaped.  If  there  is  a  favorable  pressure  gradient,  then  the -wake  diameter 
will  increase  in  the  downstream  direction  due  to  the  decreasing  gas  den¬ 
sity.  Here,  the  momentum  defect, /qu^G  j ,  is  constant  with  axial  position, 
and  A  increases  as  /qu^  decreases. 

For  such  a  cylindrical  wake,  the  solution  to  the  energy  equation  (Pr  =  1) 
is  similar  to  that  given  in  Item  1  and  can  be  written  as 
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where  Jvi3  the  enthalpy  of  the  gas  at  a  temperature  equal  to  that  of  the 
tow  cable  or  riser  line.  With  the  velocity  and  enthalpy  profiles  known, 
the  wake  characteristics  can  be  completely  determined.  The  transformed 
and  physical  coordinates  then  can  be  related  using  Equation  82  for  £  =  24 
with  the  result  that 


When  the  above  equation  is  applied  at  q  =  1,  the  actual  wake  radius  or 
thickness  can  be  related  to  the  transformed  thickness  such  that 


H.  h 

0.72  +  0.08  -t~  +  0.2C 

hi  hi 


.With  the  use  of  Equation  54,  it  can  be  further  shown  that 


(p2  =  7.56  ^1  +  0.  Ill  +  0. 


278 


h  \ 

0 

2 

fi  u 

|_JL 

S  8 

h. 

1  D 

2  * 

(97) 


(98) 


With  Equations  93  through  98,  the  solution  is  complete  with  the  proper¬ 
ties  being  clearly  specified  in  terms  of  the  forebody  boundary  layer 
characteristics  and  the  inviscid  properties  at  the  edge  of  the  viscous 
wake.  In  this  approximate  solution,  the  variation  in  properties  with  X 
is  controlled  by  the  change  in  the  inviscid  flow  field  with  X  since  the  so¬ 
lution  is  effectively  a  constant  £  solution,  and  there  is  no  meaningful  re¬ 
lationship  between  £apd  X  (except  that  £is  constant). 

4.  TURBULENT  WAKE 

In  Item  1,  the  soiuticn  for  the  viscous  wake  profiles  was  obtained  except 
for  relating  £,  the  form  parameter,  to  X.  Starting  with  the  equation  for 
the  turbulent  viscosity  coefficient 


(99) 


as  given  in  Section  III,  Item  3c,  and  where  m,  n,  and  p  arc  defined  for 
the  various  viscosity  models  (see  Table  VIII). 
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TABuE  VIII  -  RELATIONSHIP  OF 
VARIOUS  VTSCOSITY  MODELS 


Substituting  Equation  99  into  Equation  63  results  in 


Introducing  the  relationship  between  Uq/uj  and  4"  leads  to 


When  A  is  replaced  through  the  relations  hip  between  A  and  8,  it  follows  that 


After  rearrangement,  the  integrable  form  is“  finally  obtained: 


—  W! Ca*  •'e*?n**p*W **.  /  *,«**«*_• 
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Integration  is  particularly  simple  for  turbulent  viscosity  coefficient 
Model  1  where  m  =  n  =  p  =  0.  The  above  then  becomes 


dX  =  - 


tn  rr~F 

<  \'T?5  -  sm 


(101) 


Here,  Xq  is  the  point  where  the  turbulent  wake  begins,  and  £q  is  the 
value  of  i;at  that  point.  If  the  wake  is  initially  laminar,  then  Xq  is  the 
physical  location  of  the  transition  point,  and  Cq  equals  Laminar  at  tran¬ 
sition.  Integrating  Equation  101  results  in 


f  <*  -  V 


/  80  6.  5  \ 

ILL  JL 

W2'  *  / 

V 120  "  5500 

(102) 


Since  0  was  considered  constant  in  Equa'  ons  101  and  102,  then  these 
equations  are  limited  to  regions  of  zero- pressure  gradient.  Equation  101 
is  not  so  limited. 

Another  case  of  interest  is  when  Model  4  is  used  for  the  turbulent  vis¬ 
cosity  coeffivicnt.  This  is  the  same  model  used  by  Less  and  Hromas 
(Reference  50).  In  this  case,  n  -  D  =  0  but  m  =  1,  and  it  is  necessary 
to  evaluate  p\l  Pq  before^  and  X  can  be  related.  The  value  of  Pq  is  given 
by  Equation  85.  This  equation  can  be  rewritten  as 

Pn 


=  -V 

\  -  1 


(103) 


tH- 


„4\?- 

,2  7— r  T,  24  2  1  \  zij 


Thus, 


When  this  value  is  substituted  in  Equation  100,  there  results  for  this  case, 
i.  e.,  using  **odel  4 
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|dx  = 


LM. 

_±  j± 

120  ‘  5500 


y 

1  +  X-„i  M2  -  (l  +  ■■■■■  M2) - -  —  -  ^  ~  M2  (l  -  2  — 

1  2  1  V  2  1/  2  24  2  l\  24/  |£ 


(105) 


For  an  adiabatic  forebody  wall,  0Tj  equals  zero;  upon  integration,  the 
a1  ■  ■  equation  becomes 


i  <x  -  V  1  - 


(‘  '  M  |  ■  v  --1-  ul  t-t  ^  M2  illUf  (106) 

r  - -  1  Z  MI\12  2  157^2’  1 


t  jfi. 

JTZO  '  5500 


where  again  0  has  been  assumed  constant,  and  Xq  and  £q  are  as  pre¬ 
viously  defined  for  turbulent  flow. 

g2  2 

Another  limiting  condition  is  that  of  a  cold  wall  where  Tj  -  0j  .  Using 
this  assumption.  Equation  100  takes  a  final  form  upon  integration  of 


t  <*-  V  =  ' 


i  a 

y  1 2  < 


(il4L.fi  -  i-  +  m2  (i-  -  ii)W , 

- 21  2  1  \24  576/]^2 


i  t  120  5500 

C0 


where  0  also  has  been  considered  a  constant. 


(107) 


For  Model  3  and  £  small,  the  relationship  between  X  and  ^(Equation  75) 
takes  the  limiting  form 


Since  for  £  small 


X  cr  <\  . 


2  - 1 

A  <*•£ 


6  *  A, 
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then  it  can  be  shown  that 


1 

x  T 
i «  X  - 


(108) 


This  is  in  the  limit  of  small  £  or  large  X,  i.  e..  the  far  wake. 

With  Models  2  and  4,  it  can  be  shown  that  for  £  small  and  *>q  *  ^ 

dX(*  , 

Then, 


X»£Z  . 

2a  a  similar  manner,  the  limiting  farm  of 


6  *X" 


(109) 


is  arrived  at  far  the  far  wake  region.  Furthermore,  for  Models  2  through 
4,  it  also  can  be  shown  that  for  f  small 


(110) 


(HI) 


These  also  are  in  agreement  with  other  investigations  of  the  far  wake. 

Only  in  the  case  of  Model  1  is  a  deviation  noted  from  the  proper  asympto¬ 
tic  behavior  as  indicated  by  experimental  data.  With  Model  1, 

dXo<^|, 

and  thus 


X«£2. 
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<U2) 


(113) 


Thi*  lack  of  a  Jroper  behavior  in  the  limit  of  large  X  casts  some  doubt  on 
the  suitability  of  Model  1. 


(Reverse  is  blank) 
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APPENDIX  V 

WIND-TUNNEL  TESTING  AND  DATA  REDUCTION 


1.  GENERAL 

AH  wind-tunnel  tests  in  support  of  the  EUREKA  program  were  conducted 
in  the  40-in.  supersonic  tunnel  [Gas  Dynamics  Wind  Tunnel,  Supersonic 
(A)j  of  the  von  Karman  Gas  Dynamics  Facility  (VKF)  at  the  Arnold  En¬ 
gineering  Development  Center  (AEDC).  The  complete  series  is  as  fol¬ 
lows; 

1.  WT-1  -  static  stability  tests  in  support  of  design 
feasibility  study  described  in  Volume  I 

2-  WT-IIA  -  wake  survey  and  forebody  surface  pres¬ 
sure  tcet  (axi symmetric  and  asymmetric  forebodies) 

3.  WT-DB  -  parachute  performance  tests  (axisym- 
metric  and  asymmetric  forebodies) 

4.  WT-ni  -  shock  shape  determination  (axisymmetric 
and  asymmetric  forebodies);  surface  pressure  dis¬ 
tribution  (asymmetric  forcbody) 

5.  WT-IV  -  wake  centerline  survey  (axisymmetric 
forebody) 

6.  WT-V  -  wake  survey  (axisymmetric  forebody) 

7.  WT-VI  -  parachute  performance  tests  (axisymmet¬ 
ric  and  asymmetric  forebodies) 

The  tabulated  wake  survey  and  forebody  surface  pressure  data  are  given 
in  Volume  IIL  The  static  stability  tests  are  described  in  Volume  I;  the 
remaining  tests  are  described  In  tins  .ppeodix. 

t.  WIND-TUNNEL  DESCRIPTION 

Tunnel  A  is  a  continuous,  closed-circuit,  variable  density  wind  tunnel 
with  an  automatically  driven,  flexible  plate-type  nozzle  and  a  40-  by 
40-in.  test  section.  The  tunnel  operates  from  Mach  1.  5  to  6  at  maxi¬ 
mum  stagnation  pressures  from  29  to  200  psia,  respectively,  and  stagna¬ 
tion  temperatures  up  to  300  6).  Minimum  operating  pressures 

*re  about  one-tenth  the  maximum  at  each  Mach  number. 

3.  MODELS  AND  SUPPORT  SYSTEM 

Figure-  29  shows  the  forebody  models  in  detail.  The  basic  Arapaho  C  and 
the  Arapaho  C  with  nose  cone  forebodies  are  axisymmetric.  The  blunted 
elliptical  cone  forebody  has  an  elliptical  base. 
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For  the  woke  survey  aud  parachute  performance  tests,  the  forebodies 
were  supported  with  *  horizontal  wall-to-wall  strut  having  a  constant 
chord  and  a  constant  thickness.  For  the  parachute  tests,  the  riser  line 
passed  through  the  forebody  and  mounting  strut  and  into  a  plenum  tank 
where  it  was  attached  to  the  tensiometer  and  winch  assembly.  Figure  44 
shows  this  setup. 

The  models  were  sting-mounted  from  the  rear  for  the  shock  rh~pe  deter¬ 
mination  (WT-m). 

4.  INSTRUMENTATION 

From  calibration  results,  t}ie  maximum  variation  of  tunnel  centerline 
Mach  number  is  about  *0. 5  percent. 

The  tunnel  stilling  chamber  temperature  is  accurate  to  ±3  R.  The 
tunnel  stilling  pressure,  po,  was  measured  with  transducers  of  5,  15, 

30,  40,  150,  and  300  psid  capacity  that  arc  considered  accurate  to  within 
0.3  percent  of  full-scale  capacity. 

The  model  base  pressures  were  measured  with  1,  5,  and  15  psid  capacity 
transducers,  referenced  to  a  near  vacuum,  which  also  are  considered 
accurate  to  within  0. 3  percent  of  the  transducer  capacity.  The  lowest 

"  base  pressures  were  encountered  at  the  M^  =  5  low  Reynolds  number 
test  conditions.  At  the»e  test  conditions,  base  pressures  were  about 
0.008  psia.  The  0.3  percent  accuracy  of  a  1  psid  transducer  then  could 
introduce  errors  of  approximately  *37.  5  percent  of  the  measured  values. 

The  local  wake  pitot  and  static  pressures  were  measured  with  transducers 
of  the  same  capacity  ns  those  used  Tor  base  pressure  measurement.  The 
lo  est  pressures  and,  therefore,  the  greatest  inaccuracies  were  en¬ 
countered  at  th-  Mflj,  "  5  minimum  pressure  test  conditions.  Combined 
errors  in  the  local  static  and  pitot  pressures  may  introduce  inaccuracies 
of  0.2  M  at  Mj^  =  1.0.  The  inaccuracy,  in  this  case,  was  determined 
using  maximum  errors  in  pressure  measurements  assuming  (1)  the  local 
static  pressure  was  equal  to  free-streaxn  static  pressure  and  was  mea¬ 
sured  with  a  1  psid  transducer  and  (2)  the  pitot  pressure  was  measured 
with  a  5  psid  irausoocer.  Making  these  same  assumptions  tt  free- stream 
Mach  number  of  4,  3,  and  2  and  at  the  minimum  pressure  condition,  in¬ 
accuracies  in  measuring  a  local  Mach  number  of  1.0  are  0.  11,  0.  06,  and 
0.03  M^,  respectively.  At  the  high-pressure  conditions  for  these  tests, 
the  accuracy  of  measurement  at  =  5,  =  l.C  is  increased  by  ap¬ 

proximately  four.  The  lower  Mach  number  test  cases  reflect  a  similar 
accuracy  increase  at  the  high- pressure  test  conditions. 

Table  IX  lists  the  test  conditions  for  the  wake  survey  tests.  Based  on 
this  table,  the  error  in  the  measurement  of  the  local  static  pressures, 
which  are  approximately  equal  to  the  free-stream  static  pressure  and 
therefore  are  the  mpst  critical  for  measurement,  due  to  the  inaccuracy 
of  a  1  psid  transducer  appears  in  the  third  significant  digit  except  for 
the  Mg,  a  5  case.  The  total  pressures,  being  higher  than  the  static 
pressures,  then  should  reflect  a  higher  degree  of  accuracy. 
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TABLE  DC  -  WAKE  SURVEY  TEST  CONDITIONS  SUMMARY 


Re_ 

- - - i 

i 

r - 

_ 22.  X 

in.  * 

P 

P 

T 

rrODe  Itatioas 

Forebodies 

Mcc 

106 

0 

i>iria 

00 

(psia) 

o 

fdeg  R) 

X> 

S' 

Basic  Arapaho 

*» 

%> 

0.  05 

2,8 

0.360 

581 

1.  2,  3.  5,  7, 

C,  Arapaho  C 
with  nose  cona, 

2 

0.365 

19.0 

2.430 

582 

i  1 

9,  and  1 1  for 
all  teat  points 

and  blunted 

3 

0.  06 

5.9 

0.  17C  | 

I  583 

elliptical  cone* 

3 

0.  50 

43.0  | 

1.  170 

582 

f 

■: 

4 

0.09 

13.  5 

0.  090 

583 

| 

j 

i 

4 

0.  50 

73.0 

0.460 

583 

• 

'  ’j 

5 

0.13 

| 

30.  5 

0.058 

585 

1 

5 

0.  50 

133.0 

0.251 

1  631 

i 

Arapaho  C  with 

2 

0.06 

2.8 

0.  360 

560 

0.5.  1,  2,  3, 

nose  cone  only* 

3 

0.075 

5,9 

0.  170  I 

560 

5,  6.  7,  9, 
and  11  for  ell 

4 

0. 10 

13.5 

0.090 

560 

test  points 

5 

0. 135 

30.  5  ! 

0.058  i 

560  ] 

Arapaho  C  with 

■*» 

I  * 

0,06 

2.8  j 

0.360  ! 

560 

0.5,  1,  2,  3, 

nose  core  only* 

3 

0.075 

5.9  j 

0.  170 

560 

5,  6,  7,  9, 
and  1 1  for  all 

4 

0.  10 

13.  5 

0.090 

560 

test  points 

5 

!  0.135  j 

30.  5 

C.  058  j 

5b  C 

1 

WT-IIA,  pitot  and  static  rakes. 

WT-rV,  pifot,  static,  and  total  temperature  centerline  probes. 
*WT-V,  pitot  rake. 


The  wake  survey  rakes  used  in  WT-IIA  and  WT-V  were  of  a  cruciiorm 
design  with  multiple  probes,  the  reby  making  it  possible  to  obtain  complete 
wake  profiles  at  each  X/D  station  without  lateral  or  vertical  movement  of 
the  rake.  Table  X  lists  location;-  of  the  probes;  Figu-’  45  sho>vi  the 
probe  geometry. 

The  centerline  probes  used  in  WT-IV  have  the  same  geometry  as  those 
used  on  the  multiprobe  -akee. 

The  WT-ilB  parachute  urag  measurements  were  made  witn  a  200  ten¬ 
siometer  located  in  the  winen  plenum  tank.  A  ti"  nistory  of  the  dmrmc 
drag  output  from  the  tension- eter  was  recor-’.a  on  an  osciliogra ph,  and 
the  average  drag  values  were  determined  from  i’ne  recorded  traces.  The 
accuracy  of  the  200-lb  tensiometer  determined  from  static  load*  is  as 
follows. 


12C 


lKE  survey  probe  locations 


j 

i'  itcc 

WT 

-HA 

WT- VI 

ra  kr 

!  Static 

rake 

Pitot 

rake 

nu 'Tiber 

*  (in.  ) 

z  (in.  ) 

j  Y  (in.  } 

Z  (in. ) 

Y  (in. } 

Z  (in. ) 

1 

i  -1,610 

0 

|  -0.966 

0  " 

-4.  57 

c 

2 

i  0 

0 

(g!£|  v 

0 

-4.  1.7 

0 

3 

0.322 

0 

j  0.322 

0 

-3.77 

0 

0.  644 

:  0 

j  0.805 

0 

-3.37 

0 

j  0. 966 

0 

0 

-1.449 

-2.97 

0 

6 

1.  288 

0 

i 

0 

-0.966 

-2.  57 

0 

7 

!  1.610 

0 

0 

-0.483 

-2.  17 

0 

3 

1.932 

i 

0 

0 

0.322 

-1.77 

0 

Q 

1  2.  576 

0 

0 

0.805 

-1.37 

0 

10 

0 

-2. 576 

0 

a.  289 

-0.97 

0 

11 

!  o 

-1.610 

0 

1.  932 

-0.64 

0 

12 

0 

0.322 

0 

2  .  576 

-0.32 

0 

13 

0 

0.644 

2.  17 

0 

14 

0 

0.  966 

3.37 

0 

15 

0 

1.288 

0 

-4.  57 

16 

0 

1.610 

; 

0 

-4.  17 

17 

0 

1.932 

0 

-3.77 

18 

0 

2.254 

0 

-3.37 

19 

0 

2.  576 

0 

-2.97 

20 

0 

2.  898 

0 

-2.  57 

21 

0 

3.220 

0 

-2.  17 

22 

0 

3.864 

0 

-1.77 

23 

0 

4,  508 

0 

-1.37 

24 

0 

5.  152 

0 

0 

25 

0 

26 

0 

->** 

0 

28 

1 

0 

29 

c 

30 

0 

1.92 

31 

.  .  -  f 

0 

2 .  57 

Tl 

-0. 966 

0 

0,  966 

9 

T2 

0 

-0.  96#; 

0 

-0. 966 
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joad  range  (pound s) 


0  to  5 


6  to  15 
16  to  30 
31  tc  60 


Error  (pounds/ 

±0. 536 
*0.493 
it.  290 
*1.820 


These  accuracy  values  were  for  the  tensiometer  alone  and,  therefore,  do 
not  reflect  any  further  inaccuracies  caused  by  friction  in  the  chauge-of- 
direction  pulleys  in  the  system. 

The  WT-VI  drag  measurements  were  made  with  a  60-lb  tensiometer. 

The  tensiometer  was  calibrated  for  both  30-  and  60-ib  load  ranges, 

Tb-;  change- of-direction  pulleys  *&  the  system  were  equipped  with  needle 
bearings  in  place  of  bronze  bushings  as  used  in  previous  testa.  Average 
drag  values  were  monitored  with  a  low  response  servopotentiometer , 
and  a  time  history  of  the  dynamic  output  from  the  tensiometer  was  re¬ 
corded  on  an  oscillograph,  I‘he  accuracy  of  the  complete  drag  measuring 
system,  determined  from  repeated  static  loads  applied  before  the  test, 
is  shown  below. 


Load  range  (poc  ■dsj 


6  to  15 
it  1.0  30 
"  l  to  60 


±0. 15 
±0.  25 
*0.33 
±0.84 


Pa i  ..chute  perfo?  juiance  for  both  tests  was  monitored  on  two  high-speed, 
1.6 -mm  mot  o.i  picture  cameras  -  one  for  conventional  motion  pictures 
and  one  for  scnKeren  photography.  Still  photographs,  both  conventional 
and  schlieri-n,  also  were  obtained. 

5.  PARACHUTE  MODELS 

PARASONIC  parachutes,  which  have  a  constructed  shape  resembling  the 
inflated  shaped  of  a  Hyperflo  parachute  (Reference  Figure  28),  were 
tested.  The  Hyperflo  parachutes  are  constructed  with  a  Hat,  porous 
roof  and  a  truncated,  conical  nonporous  skirt.  The  construction  details 
of  the  PARASONIC  models  are  given  in  Table  XI. 

6.  DATA  REDUCTION 

a.  PARASONIC  Performance  Data 

The  data  obtained  from  the  PARASONIC  performance  te3ts  are  given  in 
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Table*  XQ  ^.nd  XIII.  In  addition  to  the  average  drag  values,  it  was  also 
desired  to  Have  average  ruaximum  end  minimum  drag  values.  These 
value*  were  obtained  from  the  recorded  oscillograph*  as  shown  below. 


Due  to  improper  oscillograph  calibration,  average  maximum  and  mini¬ 
mum  drag  values  could  net  be  determined  for  the  WT-VI  tests,  a  was 
possible,  however,  to  obtain  the  drag  oscillation  frequency. 

The  diametric  and  angular  measurements  given  in  Tables  XII  and  XIU 
ware  obtained  from  the  high-speed  motion  pictures. 

b.  Wake  Survey  Tests 

«i»Ov<nn>ia— — m— mmtmmmt am — — 

During  the  reduction  of  the  wake  survey  data  acquired  under  WT-LA,  a 
significant  amount  of  data  scatter  was  observed  in  formulating  velocity 
profiles  required  for  computation  of  experimental  wake  momentum  de¬ 
fects  (see  Section  HI,  Item  3a),  The  scatter  was  sufficient  to  came 
considerable  difficulty  in  evaluating  cxper.mei.tall.  achieved  wake  mo¬ 
mentum  defects  obtained  by  integration  of  the  velocity  profiles.  To  com¬ 
bat  this  difficulty,  the  wake  data  for  each  test  condition  were  smoothed  by 
means  of  a  multiple  regression  technique  to  permit  more  meaningful  in¬ 
tegration  of  the  velocity  profiles.  Tins  achieves  by  expressing  the 
family  of  curves  formed  by  the  velocity  profile  data  points  as  a  two-di- 
meuslpn&l  power  series  (i.  e. ,  a  pec/er  series  in  Y/D  -  or  Z/D  -  and 
their  cross  products)  to  produc-  a  best  r'-aiisilctl  fit  for  the  family  cf 
curves.  The  best  fit  use*  an  equation  of  the  form 
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TABLE  XI'  Mf  *i>EE  DECEL,e.RATOR  FAB 


rameters 


D 

open 

(in. ) 


rm 


(iru)  J  (in.)  NS  j  NG 


7.87  |  6.0 


3.06 
3.  04 
3.04 
3.05 
3.02 


6.70  j  6.0 
9.28  |  6.0 


7.86  j  4.06 
7. 86  \  4.0 

7.86  4.01 


8.90  7.86  3.70 

9.60  7.86  3.82 

8.27  7.86  3.54 


8.04 

7.64 

7.50 


7.  55 
7.50 
7.55 


12.69  12  j  12 

12.50  12  |  12 


'12.38  12 


Roof  me* 


I 


0.25  I  1.625 


2.300 

2.300 

2.300 


Materials,  coatings. 

specifications 

Roof  » *d£sh 

i 

Suspension  lines 

Roof  cap, 
skirt  panels 

Skirt  reinforcement  web, 
apex  reinforcement  web, 
roof/skirt  seam,  radial s 

HT- 83-44  Nomex  cloth 
(172-Ib  tensile,  warp) 
3.01  ©*/»q  Pre- 
coat  ed  with  D1569-. 

F839  polyurethane, 
porosity  controlled  by 
coating  with  0-65  flex¬ 
ible  thermal  coating 

M1L-C-5Q40  Type  I 
cord,  100-  lb  tensile 
strength 

0. 84  oz/ sq  yd 
4787  nylon 
coated  with 
U1569-F839 
polyurethane 

Double  fold,  cotton  bias 
tape,  88  X  80  count, 

20  Ib/in. 

MUL-T-713A,  Class  H. 
Type  8  nylon  lacing 
(50- lb  tensile)  pre¬ 
coated  on  form  with 
D1569-F839  polyure¬ 
thane,  porosity  con¬ 
trolled  by  coating  with 
D-65  flexible  thermal 
coating 

MIL-C-5040  Type  I 
cord,  100- lb  tensile 
strength 

0.84  oz/sq  yd 
4787  nylon 
coated  with 
D1569-F839 
polyurethane 

Double  fold,  cotton  bias 
cape,  88  X  80  count, 

20  Ib/in. 

HT-83-44  Nomex  mesh 
cloth  (172-lb  tensile, 
warp)  3.01  oz/sq  yd, 
precoated  with  D1560- 
F839  polyurethane, 
porosity  controlled  by 
coating  with  D-65  flex¬ 
ible  thermal  co  ting 

MIL- C- 5040  Type  m 
cord  with  fibers  re¬ 
moved 

0. 84  oz/sq  yd 
4787  nylon 
coated  with 
D1569-F839 
pol /urethane 

Double  fold,  cotton  bias 
tape,  88  X  80  count, 

20  lb/in. 

MIL- C- 5040  Type  I 
wiui  fibers  removed 
precoated  on  form 
with  D1569-F839 
polyurethane,  porosity 
controlled  by  coating 
with  D-65  flexible 
thermal  coating 

MIL-C-5040  Type  HI 
cord  with  fibers  re¬ 
moved 

0.  84  oz/ sq  yd 
4787  nylon 
coated  with 
D1569-F839 
polyurethane 

Double  fold,  cotton  bias 
tape,  88  X  80  count, 

20  Ib/in. 

(Reverse  is  blank) 
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Y/D  Case* 


Z/D  Cases 


The  curve  fits  used  only  even  powers  of  Y/D  or  Z/D  because  of  symmetry 
tha*  exists  in  the  wak  '  of  the  forebedy. 

The  curve  fits  are  made  to  the  experimentally  derived  velocities  and  are 
the  best  statistical  fits  to  all  the  data  presented.  Thus,  the  curve  is 
matched  to  the  data  points  with  scatter  as  well  as  to  the  more  accurate 
da**  points.  The  multiple  regression  solution  then  gives  the  least  squares 
value  of  the  coefficient#  for  a  particular  sample  of  observations. 

Samples  of  the  curve  fits  are  shown  in  Figures  46  through  50.  The  si¬ 
nusoidal  type  variation  of  the  polynomials  above  Y/D  values  of  0.4  are 
smoothed  out  by  the  integration  since  the  area  under  the  curve  is  approxi¬ 
mately  the  same  as  would  exist  without  the  sinusoidal  type  variation  of 
velocity  predicted  by  the  polynomial.  The  sinusoidal  trace  is  caused  by 
the  polynomial  curve  fit  as  the  polynomial  tries  to  pass  through  the 
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Figure  48  -  £xperi'  atttal  Wake  Velocity  Data  Curve  Fit,  X/D  *  7 
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"*ig«ra  50  -  Experiment*^  ISTfck .*  Velocity  Dti*  Curve  Fit,  X/D  *  11 
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existing  data  points.  Where  the  velocity  is  constant  or  e  pproximately 
constant  with  the  change  in  Y/D  or  Z/D,  the  polynomial  is  unsuccessful 
in  matching  the  experimental  data.  The  one  polynomial  for  each  case  i. 
the  oest  lit  to  all  the  data  points  at  each  X/D  stauon;  therefore,  the  poly¬ 
nomial  is  attempting  to  fit  or  match  a  surface  of  points  in  two  dimensions 

c.  Test  Data 


Tables  XII  and  XIII  give  the  parachute  performance  test  data  (Table  XII 
has  data  lor  tne  WT-IIB  seri  s  and  Table  XIII,  for  the*  WT-VI  scries). 
Ectb  dr  a"  and  li^v  nsrameten  weu  determined. 


(Reverse  ;•  blann) 
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Methods  were  developed  to  predict  the  boundary  layer  and  wake  characteristics  of 
symmetric  bodies  moving  at  s'lpeescnic  speeds.  The  methods  developed  for  predict¬ 
ing  boundary  layer  characteristics  are  applicable  to  both  laminar  and  turbulent  flow 
over  axisymmetric  or  two-dimensional  bodies  and  can  be  used  as  an  approximation 
for  the  flow  over  quasisymmetric  bodies.  The  wake  predictive  methods  were  de¬ 
veloped  for  axisymmetric  fiow  and  are  applicable  to  both  laminar  and  turbulent  wakes 
The  aforementioned  methods  were  combined  and  programmed  so  that  wake  profiles 
can  be  determined  based  only  on  inputs  of  forebody  geometry  and  free  stream  con¬ 
ditions. 


Also  presented  is  a  method  for  calculating  the  drag  of  PARASCNIC  parachutes  using 
the  flow  field  immediately  ahead  of  the  parachute  as  the  free -stream  conditions.  An 
analysis  of  the  effects  of  various  geometric  and  free -stream  parameters  cn  PARA- 
SONIC  parachute  performance  was  conducted  and  trends  determined. 
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